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Abstract. We prove explicit formulas for the p-adic L- functions of totally real number fields 
and show how these formulas can be used to compute values and representations of p-adic 
//-functions. 

1. Introduction 

The aim of this article is to present a general method for computing values and represen- 
tations of p-adic L-functions of totally real number fieldsQ These functions are the p-adic 
analogues of the "classical" complex L-functions and are related to those by the fact that they 
agree, once the Euler factors at p have been removed from the complex L-functions, at nega- 
tive integers in some suitable congruence classes. The existence of p-adic L-functions was first 
established in 1964 by Kubota-Leopoldt [22J over Q and consequently over abelian extensions 
of Q. It was proved in full generality, 15 years later, by Deligne-Ribet [T2] and, independently, 
by Barsky [3] and Cassou-Nogues [7J. The interested reader can find a summary of the history 
of their discovery in [7j. 

There have already appeared many works on the computation of p-adic L-functions, starting 
with Iwasawa-Sims |20| in 1965 (although they are not explicitly mentioned in the paper) to 
the more recent computational study of their zeroes by Ernvall-Metsankyla |16[ [T7] in the 
mid-1990's and the current work of Ellenberg-Jain-Venkatesh [15] that provides a conjectural 
model for the behavior of the A-invariant of p-adic L-functions in terms of properties of p-adic 
random matrices. However, most of these articles deal only with L-functions over Q or that 
can be written as a product of such L-functions. One remarkable exception is the work of 
Cartier-Roy [6] in 1972 where computations were carried on to support the existence (at the 
time not yet proven) of p-adic L-functions over some non-abelian fields of degree 3, 4 and 5. 

The method for computing p-adic L-functions given in the present paper is derived from the 
construction found in [JJ [21] [23] . It generalizes a previous work with Solomon [27] . The idea 
is the following. First, using Shitani's cone decomposition (see Subsection I3.3p . we express 
L-functions in terms of cone zeta functions (see Subsection 13.41 Proposition 13.11 and Equa- 
tion [3]6]). Then, for a given cone zeta function, its values at negative integers are encoded 
into a power series (see Subsections 13.41 and I3.5p . Using the method of Section [2j this power 
series is then interpreted as a p-adic measure. The p-adic cone zeta function is obtained by 
integrating suitable p-adic continuous functions against this measure (see Theorem I2.9P once 
it is proved that it satisfies the required properties (see Subsection 14. ip . The main tool for the 
computation is an explicit formula for the power series associated with the cone zeta function, 
up to a given precision; see Theorem 14.11 From this formula we give an explicit expression 
for the values of the cone zeta function at some p-adic integer (Theorem 15. 28p and an explicit 
expression for the corresponding Iwasawa power series (Theorem I5.24p . Note that these also 
are valid only up to a given precision. 
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One shortcoming of our method is that it is not very efficient compared to the complex 
case (see Subsection 15.51 for some complexity estimates). For example, in this simplest case 
of p-adic L-functions over Q, for a Dirichlet character of conductor /, the complexity in / 
of the method presented here is 0(/ 1+e ), whereas there exist methods to compute complex 
Dirichlet L-functions in 0(f 1 ^ 2+t ). Even in this simple case it remains an open problem 
whether methods as efficient exist in the p-adic case. 

Note. The construction presented in this paper was developed over several years and during 
that time was used in two previous works; see [28] . It is worth noting that the method has 
evolved and therefore the brief description of it in these earlier articles does not necessarily 
match exactly the one that is finally presented here. 

Acknowledgments. I would like to express my gratitude to D. Solomon for his invaluable 
help with this project and for introducing me to the fascinating world of p-adic L-functions. 
I am also grateful to D. Barsky for his help at the start of the project and to R. de Jeu and 
A. Weiss for providing with the motivation to complete this project. Finally, I thank heartily 
D. Ford for reading carefully the earlier version of this article and for his many corrections 
and helpful comments. 

2. p-ADIC INTERPOLATION 

Let p be a prime number. Denote by Q p the field of rational p-adic numbers. The subring 
of p-adic integers is denoted by Z p , and C p is the completion of the algebraic closure of Q p . 
Let | • | p denote the p-adic absolute value of C p normalized so that \p\ p = p~ l and v p (.) the 
corresponding valuation; thus v p (p) = 1. For / > 1, an integer, let Wf denote the subgroup 
of f-th roots of unity in C p . The torsion part T p of the group Z* of units in Z p is equal to 
WtpU\ where q := 4 if p = 2, q := p if p is odd, and tp is Euler totient function. We have 
Z* = T p x (1 + qZ p ), and the projections to : Z* — > T p and (•) : Z* — > 1 + qZ* are such that 
x = uj(x){x) for all x £ Z* . In particular, we have x = uj(x) (mod q) for all x G 



'X 



2.1. Continuous p-adic functions. For n G N := {0, 1, 2, . . . }, the binomial polynomial is 
defined by 

'l ifra = 0, 

x(x-l)---(x-(n-l)) , (2.1) 
n J — ^ i otherwise. 

The binomial polynomial takes integral values on Z, hence, by continuity, it takes p-adic 
integral values on Z p . Let / be a function on Z p with values in C p . One can easily construct 
by induction a sequence (f n )n>o of elements of C p (see Subsection 15. 1|) such that 

f(x) = Yjn( X ) forallxGN. (2.2) 



n>0 v 7 

(Note that this is in fact a finite sum.) The coefficients f^s are uniquely defined and are called 
the Mahler coefficients of /. We have the following fundamental result (see \26\ §4.2.4]). 

Theorem 2.1 (Mahler expansion). Let f be a function on Z p with values in C p . Then f is 
continuous on 7L p if and only if 

lim \f n \ p = 0. 

n— >oo 

If f is continuous, then the sequence of continuous functions 

N 



*h>5>(;) (2.3) 

n=0 V ' 

converges uniformly to f . Reciprocally, let (f n )n>o be a sequence of elements inC p converging 
to zero. Then the sequence of functions in ()2.3[) above converges to a continuous function. 
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Denote by C(Z p ,C p ) the set of continuous functions on Z p with values in C p . For f G 
C(Z p ,C p ), we define the norm of / by 



P ■= max|/(x)| p . 

The norm of / is a finite quantity since Z p is compact and in fact, if (f n )n>o are the Mahler 
coefficients of /, we have 

||/|| p = max|/ n | p . (2.4) 

n>0 

2.2. A family of continuous functions. We define a family of continuous functions that 
will be useful later on. For s G Z p , we would like to define x \-t x s , where x is an p-adic 
number, in such a way to extend the definition of x \— > x k when s = k G Z. In general, it is 
not possible. However, when x G 1 + </Z p , one can set 



n>0 

The series converges since |x — l| p < 1, and by Theorem l2.ll the function s i— )• x s is continuous^ 
Furthermore, when s = G N we recover the usual definition of x fc by the binomial theorem, 
and it is easy to see that this function has the expected properties. With that in mind, for 
s G Z p we define the function <p s in C(Z p , C p ) by 



<f> s (x) :- 




if x G pZ p , 
if x G Zp . 



It it easy to see that <j) s is a continuous function and that its restriction to Z^ is a group 
homomorphism. We state two results concerning the properties of <p s . The first one follows 
directly from construction. 

Lemma 2.2. Let k be a integer. Then, for all x G 1 + gZ p , we have 

4> k {x) = x k . □ 
Lemma 2.3. The map /rom Z p to C(7j p ,C p ) is continuous. 

Proof. Let s, s' in Z p . If x G pZ p , then (f) s (x) = (j) s '(x) =0. If x G Z^ , then 

|0 s (x) - y (a:)| p = |(x) s (l - (x) s '-) 

(x - 1)™ 



v- is' - s 



n>0 



(x - l) n 

n 



< Is' 



< |s — s\ p . □ 

p 

2.3. Integration of p-adic continuous functions. A measure fi on Z p is a bounded linear 
functional on the <C p -vector space C(Z P , C p ). That is, there exists a constant B > satisfying 

Hf)\ P <B\\f\\ p for aU/€C(Zp,Cp). (2.5) 

The smallest possible -B is called the norm of the measure /x and is denoted \\fJ>\\ p . With this 
norm, the set A / l(Z p ,C p ) of measures on C(Z p ,C p ) becomes a C p -Banach space. From now 
on, we will write 

J f(x)dfi(x) :=fi(f). 
Usually we will drop the x to simplify the notation when the context is clear. 
Lemma 2.4. The function fx is a continuous map from C(Z P ,C P ) to Z p . 



2 Actually, for p = 2 we only need x G 1 + 2Z2 . But in order to have an analytic function it is necessary to 
assume x G 1 + 42^ ; see Subsection 15.21 
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Proof. This is clear by (|2.5p . □ 
Lemma 2.5. Let f G C(Z p ,C p ) toit/i Mahler coefficients (f n )n>o- Then we have 

[fdfi = J2fn[ (^W (2.6) 

J n>0 ^ 

N fx\ 

Proof. This is clear since / = lim V f n \ and is continuous by the previous lemma. □ 

iV->oo n=0 \nj 

For /i G .A/f(Z p ,C p ) and 5 G C(Z p ,C p ), the measure g/i is defined, for any / G C(Z P ,C P ), by 

f(x)dgn(x) := / f(x)g(x)dfj,(x). 



When g = xa> the characteristic function of an open and closed subset A of Z p , we will use 
the notation 

fdfi:= J fdxAf-i- 

A measure \i is said to have support in A ii fj, = XAfJ-- In other words, for all / G C(Z p , C p ) we 
have 



//*=-/ 



fdfi. 



2.4. Measures and power series. Let C P [[T]] M be the C p -algebra of power series whose 
coefficients are in C p and are bounded in absolute value. Let \x be a measure in .A/f (Z p , C p ). 
One associates with \i a power series G C P [[T]] M defined by 



n>0' 

l6d 



Reciprocally, given a power series F G C P [[T]] M with coefficients F n (n > 0), one associates 
with F a measure /zp defined by 



E / ( J ■= F «- ( 2 - ? ) 

n>0*' 

Indeed, by Lemma 12.51 these equations uniquely determine the measure /if. These maps 
define isometric isomorphisms of C p -Banach space between M(Z p ,C p ) and C p [[T]] bd where 
the norm on C p [[T]] bd is defined to be the maximum of the absolute values of the coefficients; 
see [231 Chap. 4]. 

Remark 2.6. Another characterization of the correspondence between measures and bounded 
power series is that the power series F^ is the unique power series in C p [[T]] bd such that 

J (1 + t) x dfi(x) = Ffj,(t) for all t G C p such that \t\ p < 1. 

The measures corresponding to powers of 1 + T form an important class. The result below 
follows directly from (j2.3jl . Lemma 12.51 and (|2.7p (or simply the remark above). 

Lemma 2.7. Let a G Z p . Then the measure associated with the power series 

(l+T) a :=E( a V" ( 2 - 8 ) 

n>0 

is the Dirac measure at a, that is, the measure \i a such that 

fdfia = f(a) 



COMPUTING p-ADIC L-FUNCTIONS OF TOTALLY REAL NUMBER FIELDS 

for all f £ C(Z p ,C p ). 



5 

□ 



d 



2.5. The interpolation principle. Let A be the linear operator (1 + T)--^ acting on 
C p [[T]] bd . Let [i E M(Z p ,Cp) be a measure. We have 



AF,(T) = (1 + D £ n / ( X ) d„ ■ T n ~ x = W 

n>l J ^ ' n>0 



n>0 



n>l 



(n + 1) 



n+ 1 



+ n 



dfi-T n = F x JT). 



Thus we have proved the first part of the result below; the second follows from Remark 12.61 

Lemma 2.8. Let n be a measure with associated power series Fa- Then the measure associated 
with the power series AF^ is x^i. In particular, 



A fe F M (T)| T=0 = 

for any integer k > 0. 

We can now state the main result of this section. 



x k dji 



□ 



Theorem 2.9. Let (a n ) n >o be a sequence of elements of C p . Assume there exists a power 
series F G C P [[T]] M such that for all k > we have 

A k F(T) lT=0 = a k 

and that the associated measure fj,jr has support in 1 + qL v . Let f : Z p — > C p be defined by 



f(s) := / (f> s (x)d/j lF (x). 



Then f is a continuous function such that 

f(k) = a k 

for all k G N. 

Proof. It is clear from Lemmas 12.31 and 12.41 that / is continuous. For k E N we compute 



(2.9) 



/(*) 



4> k {x)diJ, F (x) 



x dm 



since [i has support in 1 + qL p 
by Lemma 12.21 



since \i has support in 1 + qZ p 
by Lemma 12. 



X d(J,F 



= A F(T)\ T=0 
= ak- 

3. Values of zeta functions at negative integers 



□ 



Let E be a totally real number field of degree d with ring of integers Z#. We consider E, 
and all other number fields, as subfields of the algebraic closure Q of Q contained in C. We 
also fix once and for all an embedding of Q into C p . For a in E, we denote by a® G M, 
i = l,...,d, its conjugates. An element a G E is totally positive if oft' > for i = 0, . . . , d. 
We write a ^> 0. The subgroup of totally positive numbers in E x is denoted E + and we let 
Z^ := E + n 7Le be the set of totally positive algebraic integers in E. Let M = Me/q denote 
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the absolute norm of the group 1(E) of ideals of E. By abuse, for a a non-zero element in E 
we write N(a) := N(olLe)- When a is totally positive N(a) equals the absolute norm of a. 

Let m := fj be a modulus of E, that is, the formal product of an integral ideal f of E (the 
finite part) and a subset 3 of the set of infinite places of E (the infinite part). We use the 
notations 

E m for the subgroup of elements of E that are congruent (multiplicatively) 
to 1 modulo m, 

I m (E) for the subgroup of fractional ideals of E that are relatively prime to f, 
Cl m (E) for the ray class group of E modulo m, that is, the quotient of I m (E) by 
the subgroup of principal ideals generated by elements of E m , and 

h m (E) for the cardinality of Cl m (E). 
Finally, we set U m (E) := U(E) n E m where U(E) is the unit group of E. 

3.1. Twisted partial zeta functions. Let a be a fractional ideal of E, relatively prime to 
f. The partial zeta function is defined, for s G C with 3t(s) > 1, by 

where the sum is over all the integral ideals b that are in the same class of C\ m (E) as the 
inverse of a. For c, an ideal of E, relatively prime with f, the twisted partial zeta function is 
defined, for M(s) > 1, by 

Z m (a-\qs) :=M(c) 1 - s Z m ((ac)- 1 ;s)-Z m (a- 1 ;s). (3.1) 

The partial zeta functions have meromorphic continuation to the complex plane with a 
simple pole at s = 1. Since the partial zeta functions all have the same residue at s = 1, they 
cancel out in (|3.1|) and the twisted partial zeta functions have analytic continuation to the 
whole complex plane. 

Let x be a character on the ray class group Cl m (E), and let L m (x; s) be the corresponding 
Hecke L-function defined, for 3?(s) > 1, by 

Ln(x;s):=H(l-x(qW(<\)-T 1 ( 3 - 2 ) 

qfm 

where the product is over all the prime ideals of E not dividing the finite part f of the modulus 
m. The link between Hecke L-functions and partial zeta function provides a way to express 
the former in terms of twisted partial zeta functions. 

Proposition 3.1. For all s G C we have 

h m {E) 

( X (cW(c) 1 - s -l)L m (x;s)= £ x(a t )Z m (a-\c;s) 

i=i 

where the sum is over ideals Oj representing all the classes ofC\ m (E). 
Proof. We have 

h m (E) h m (E) h m {E) 

^ x(a i )Z m (ar 1 ,c;s)=M(c) 1 - s ]T x( ai ) Z^c)" 1 ; s) - ^ X (ck) Z^ 1 ; s ) 

i=l i=l i=l 

h m {E) 

= ( x (c)Af(c) 1 - s -l) Y, X(a i )^ m (ar 1 ;s ) 
i=i 

= ( x (c)Af(c) 1 - s -l)L m (x;s). □ 
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We now make some important additional hypotheses. 
Hypotheses. 

(HI) the finite part f of the modulus m is divisible by q; 

(H2) the infinite part 3 of the modulus m contains all the infinite (real) places of E; 
(H3) c is a prime ideal of residual degree 1. 

We will denote by c the prime number below c; therefore c = A/"(c) by (H3). 

Remark 3.2. If m does not satisfy both (HI) and (H2) we can enlarge the modulus so that it 
does satisfy these conditions and we can lift \ to a character of the new modulus. Adding all 
the infinite places to 3 to satisfy (H2) does not actually change the L-function. Replacing f by 
the 1cm of f and q to satisfy (HI) has the effect of removing the Euler factors of prime ideals 
above p in (j3.2|) . This is necessary to be able to do the p-adic interpolation. Another way to 
achieve this would be to drop (HI) and to require instead that only the elements coprime to 
q are kept in the cone decompositions (see Subsection I3.3() . From a computational point of 
view these two possibilities are basically the same. 

Remark 3.3. The construction of Cassou-Nogues [7] additionally requires c to be relatively 
prime to the co-different of E. However, as we will see in the next subsection, this is not 
actually necessary. 

3.2. Additive characters modulo c. An additive character modulo c is a group homomor- 
phism £ from the additive group TLe to the multiplicative group C x whose kernel contains c. 
We denote by X(c) the set of all these characters. 

Lemma 3.4. X(c) is a finite group of order c and all elements inX(c) but the trivial character 
have kernel c. Furthermore, for i£Zg we have 



Proof. Let x be a non-trivial character in X(c). Then Z^/Ker(x) is a quotient group of 
Z^/c = Z/cZ and thus £ is completely determined by its value on 1, which can be an arbitrary 
non-trivial c-th root of unity. Therefore there are c — 1 non-trivial characters and each has 
kernel c since the non-trivial c-th roots of unity all have order c. The last statement is the 



Proposition 3.5. Let a be an integral ideal coprime to fc. Let A be a set of representatives 
of the elements ofa(~)E m under the (multiplicative) action ofU m (E). Then for 3?(s) > 1 



Proof. An ideal b is equivalent to a modulo m if and only if there exists a £ E m such that 
b = aa~ l . Furthermore, b is integral if and only if a belongs to n E m , and two elements a 
and a' of a n E m yield the same ideal ao" 1 = a' a -1 if and only if there exists a unit e G U m 
such that a = ea! . Therefore there is a one-to-one correspondence between the integral ideals 
equivalent to a -1 modulo m and the elements of A. Thus we have 




classical orthogonality relation for characters. 



□ 




(3.3) 



Z m (a- 1 ; s) = M{aoT l )- s = M{a) s £ AA(a)" s . 



(3.4) 




We now compute 



-1 



c;s) =AA(c) 1 - s Z m ((ac)- 1 ; S 



)-Z m (a- 1 ;s 



) 
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>■ n,CL A Or A 



In ) 

a6Anc aeA 

using (|3.4p and the fact that ^4 n c is a set of representatives of ac n E m modulo U m (E), since 
a and c are coprime. Finally, we obtain the conclusion using Lemma 13.41 □ 

3.3. Cone decomposition. Let j3, Ai, . . . , X g be elements in o n Zg, with 1 < g < d, such 
that the Aj's are linearly independent. We define the discrete con4j with base point (3 and 
generators Ai, . . . , X g as the following subset of a n Z^ : 

C(P;X 1 ,...,X g ) := Ip + ^mXi with m € N for 1 < % < g \. (3.-V) 
^ i=i ' 



Following the work of Shintani [29], we have the following result of Pi. Cassou-Nogues. 

Theorem 3.6 (Cassou-Nogues). There exists a finite family {C\, . . . , C m } of disjoint discrete 
cones of E with base points belonging to aC\E m and generators belonging to (afnZg) such that 
a set of representatives of a(lE m under the action ofU m (E) is given by the union CiU - • -UC m . 

Proof. This is essentially [7, Lemma 1]. □ 

A finite set {C\, . . . , C m } of cones satisfying Theorem l3.6l is called a cone decomposition of a 
modulo m. A cone C is c-admissible if none of its generators belong to c. A cone decomposition 
{Ci, . . . ,C m } is c-admissible if all the cones Cj are c-admissible (see Remark 15.291 on the 
existence of such a decomposition). From Proposition 13.51 we have, for 5t(s) > 1, 

m 

Z m {a-\t ]S )=N{a) s E E E £(a)A/"(ar. (3.6) 

CeX(c) j=l a&Cj 

3.4. Cone zeta functions. Let C := C{f3; Ai, . . . , X g ) be a c-admissible cone and £ be a 
non-trivial element of X{t). We define the zeta function of the pair (C, £), for K(s) > 1, by 

Z(C,i-s):=Y,^)N{a)- s . (3.7) 

We associate with the same data a power series F(C,£;Ti, . . . , T^) in Q[[T11 := Q[[T\, . . . , TJ] 
in the following way. First, for r £ C, we define a power series in C[[T]] bjfl 

( 1+T ) r : =Et) T " 

n>0 V ' 

Then, for a £ E, we define the following power series in Q[[T]] 

d 

U+2T :=n( 1+T *r w - ( 3 - 



y(i) 

[1 + ±i)~ 

8=1 

And finally, we set 



G(C,£;T) := - ^ =L . (3.9) 

n(l-e(Ai)(l+T)^) 



'We will say simply cone. 

This, of course, gives the usual definition when r 6 N. 
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Remark 3.7. From Lemma 13.41 and the fact that C is c-admissible it follows that £(Aj) is a 
non-trivial c-th root of unity for all i's and thus the constant term of the denominator is non- 
zero. Therefore G(C,£;T) is indeed a power series. In fact, its constant term is an algebraic 
integer divisible only by primes above c. 

The cone zeta function of C (twisted by c) is defined by 

Z(C,c;s):= £ Z(C,Z;s), (3.10) 
£eX(c) 

and in a similar way 

G(C,c;T) := G{C,Z;T). (3.11) 

Finally, we define the A- operator acting on Q[[T]] by 

A:=II(l + 2i)^r. (3.12) 

i=l 1 

Theorem 3.8 (Shintani). The function Z(C,c;s) admits an analytic continuation to C, and, 
for any integer k > 0, we have 

Z{C, c; -k) = A h G{C, c;T)\t=o- (3.13) 
Proof. We use the following lemma from Colmez Lemma 3.2]. 

Lemma 3.9. For z\, . . . , z d E R + let f(z±, . . . , z d ) be a C°° -function such that it and all its 
derivatives tend to rapidly at infinity. For all (si, . . . , s d ) E such that Sft(sj) > for 
i = 1, . . . , d define the function 

M(f; Sl ,...,s d ):= [ f(z 1 ,...,z d ) r( Z \"'^ • • • (3.14) 

J(R+*)d T(sx) ■ ■ ■ T{s d ) z\ z d 

Then M(f; •) admits an analytic continuation to C d and satisfies 

- ( d \ ki 

M(f) -fa, -k d ) = Y[ [-g^J fi z U-- ■ > Zd)\zi=-=z d =0 

i=l ^ % ^ 

for all ,k d ) GN d . □ 

Let £ G X(c) with £ ^ 1. We define a function by 

. . . , z d ) := G(C, £; e" 21 - 1, ... , - 1) = g (3.15) 

n (1 - fMe-W) 

i=l 

where, for a G Ze and z := (zi, . . . , z d ) G C d , we set T^(a) := z\oS l "> + • • • + z d a^ . It is clear 
that this function satisfies the hypothesis of the lemma. Furthermore, for z_ G (IR + *)™ and 
a>0we have < e~ T -( a ^ < 1, and therefore, by expanding the numerator and by (|3.5p . 

ft( Zl ,. ..,z d )=Y,Z(P + n ^ + • • • + n g \ g )e- T ^ +n ^ + - +n ^ = 
Thus we find that 



, tt (l) <fel f c - Zd gW dZ d 

r(si) zi y M +* r(s d ) z d ' 
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Since 



dz 

e~ za z s — =T(s)a- s 



for a £ R + *. it follows that 



M(fc s,...,s) = J2 £(«)(a«r ■ ■ ■ (« W ) S = Z(C,b s) 
for 5R(s) > 1. Lastly, we find that 

-q^Mzu ...,z d )=((l + Ti)^fG) (C, £; e" 21 - 1, ... , e"** - 1) 
for all i's. The conclusion now follows from the lemma, ()3.10j) . and (|3. 1 1[) . □ 

Recall that we have embedded Q into C p . Thus we can see G(C, c;T) as having coefficients 
in Cp and, generalizing what we did in the first part to higher dimensions, we can try to 
interpret G(C, c;T) as a measure over Zp\ There are two problems. First, the power series 
G(C, c;T) having several variables complicates things, at least from a computational point of 
view; it would be much easier to deal with a one-variable power series^ Second, this power 
series might not have bounded coefficients, since has arbitrarily large p-adic absolute value 
when (3 is not a rational p-adic integer. This problem can be solved by making a change of 
variable in G(C, t;T), as in [3D], to transform it into a power series with bounded coefficients. 
In the next subsection we will see how to transform the power series G{C, c; T) into a one- 
variable power series satisfying the direct analog of Theorem 13.81 It will turn out that this 
power series has p-adic integral coefficients, thus solving both problems at the same time. 

3.5. The ft operator. We now explain how to construct an operator that sends G(C,c;T) 
to an one- variable power series satisfying properties analogous to that of Theorem 13.81 From 
the definition of A, we see that ft should satisfy 

fi((l + Ti) ai • • • (1 + T d ) ad ) = (1 + T) ai - ad . 

Writing T°" % = ((1 + Tj) — l) ai and developing, we get the following formal definition. Let ft 
be the linear function from C[T] to C[T] defined for (at, . . . , aj) £ N d by 

ai a d 



n(T^ ■ ■ ■ T2 d ) := ■ ■ ■ ( II^ 1 )" 1 ( ^ ) ) t 1 + T ) ni ' 

ni=0 n d =0 M=l \ n «// 



■«rf 



The lemma below establishes that this application can be continuously extended to an appli- 
cation from C[[T]] to C[[T]]. 

Lemma 3.10. Let (a x , . . . ,a d ) G N '. Then, fl(T^ ■ ■ ■ T% d ) is divisible by r max ( ai >-> a <*). 
Proof. Assume without loss of generality that is the largest of the a^s. We have 



a± o,d-i , s 

ft(T{ x ■■■T^ d ) = (-l)«i+-+ a d-i ^ •••J^ (-l) 711 ^"^" 1 ( ai ) • 

— n r, . . — n \ n l/ 

a d / \ 

(_;ga d -n d I ^ j ((1 + y^i-nd-i^d 

n d =0 

til ^d— 1 / d— 1 / \ V 

(_ 1) a 1+ ...+a d _ 1 £ ... £ (n(-l) ni ( ai ((1 + ^r- 

„,_n „, ._n\,--i \ n iJJ 



ni=0 n d _i=0 

",; 

X 

n d =0 

ai id-i / a!— 1 



nd-i 



ni=0 n d _i=0 8=1 

and every term in this sum is divisible by T ad . □ 



This is not a big problem though; in fact the computations done in |27| use two- variable power series. 
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We now prove the main property of the £1 operator, that is, that it "commutes" with the 
operator A. 

Proposition 3.11. For A G C[[T]] we have 

n(A(A)) = A(n(A)). 

Proof. By linearity and the fact that the operators A, A, and f2 are linear and continuous!! 
it is enough to prove the result for monomials T" 1 • • • T? d with (a\, . . . , a^) G N d . But any 
such monomial can be written as a finite linear combination of (1 + T\) bl ■■■(! + Td) hd with 
(&i, . . . , bd) d G N, for which the result is direct by construction. □ 

Let C be a c- admissible cone. We set 

F(C,c;T) :=0(G(C,c;T)). 

Using Proposition 13.11} the next result is a direct consequence of Theorem 13.81 

Theorem 3.12. For any integer k > we have 

Z{C,c;-k) = A k F{C,c;T)\ T=0 . □ 

Let V := {Ci, . . . , C m } be a c-admissible cone decomposition of a modulo m. We define 

m 

F m (a,c;T):=Y,F(Cj,c;T). (3.16) 

Corollary 3.13. For any integer k > we have 

Z m (a-\ c; -k) = AA(a)~ fc A fc F m (a, c; T) [T=0 . (3.17) 
Proof. Clear from (|3.6|) . □ 

To conclude this subsection, we prove that the power series F m (a, c;T) does not depend on 
the choice of the cone decomposition D. Indeed, the previous corollary prescribes the values 
of A k F m (a,qT) lT=0 for all k > which, using the following result, ensures the unicity of 
F m (a,c;T). 

Lemma 3.14. Let F(T) G C[[T}} and let k > be an integer such that T k \ AF(T). Then 
T k+i | f(t) - F(0). Furthermore, if A k F(T)\ T=0 = for all k > then F(T) = 0. 

Proof. Write F(T) := £ f n T n . We compute 

n>0 

AF(T) =f 1 + J2{(n+ l)/ n+1 + nf n )T n . 
n>l 

From this it is easy to see that T k | AF(T) implies /i = • • • = fk+i = 0) which proves the first 
assertion. For the second, let k > 1. Since A k F(T)\ T=0 = 0, that is, T | A(A fc_1 F(T)), we 
have T 2 \ A k ~ l F(T) by the first part, as A k ~ 1 F(T)\ T=0 by hypothesis. Repeating this process 
(and using the fact that F(0) = A°F(T)| T=0 = 0), we eventually get T k+1 \ F(T). Since k is 
arbitrary, it follows that F(T) = 0. □ 

4. p-ADIC L-FUNCTIONS 

We put together the results of the last two sections to construct the p-adic L-functions. 
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We leave to the reader the verification that A and A are continuous. 
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4.1. Some properties of F m (a, c;T). Let a be an integral ideal coprime to c and m. We 
prove that the power series F m (a, c; T) possesses the properties required to apply Theorem 12.91 
We start by proving a useful expression for F{C, £; T) modulo powers of T. 

Theorem 4.1. For integers k and K with < k < K define the rational function 

B k:K (x) := (-l) k Y Q (^t)" € ( 41 ) 

Let C := C(/3; X\, . . . , X g ) be a c-admissible cone and let £ be a non-trivial element of X{z). 
For N>0 define the polynomial F N (C,£;T) G Q(0[ T ] by 

(N-l)d g 
F N (C, £; T) := A(C, £ (1 + T)^+^ ]J B^^M^)) 

kl,...,kg = 1=1 

where k ■ A := &1A1 + ■ ■ ■ + k g X g £ and 

m& : -5—^ — ■ 
na-«Ai)) 

i=l 

Then 

F(C,Z;T)=F N (C,t;T) (mod T N ). 
Proof. To simplify the notation we write at := £(Aj)/(l — £(Aj)) and A := A(C,£). We compute 

ai c, m - , m{l+z)B — ^— £±S? 



n a - «Ai)(i + t)a.) n a - * «i + x) A " - 1)) 

i=l i=l 

= a (1+^ £ n a r(( l +^) Ai - 1 ) ! - 

m,...,n s >0 i=l 

Let X be the ideal of C P [[T]] generated by the monomials T" 1 • • • T^ d where max(ai, . . . , cy) > 
iV. For any P(T) G X, it follows by Lemma EHD] that ft(P) G T^CpffT]]. Furthermore, for 
i = 1, ...,(? we have 

((l + T) A '-l)"*el if ni > (N - l)d + 1. 

It follows that 

(N - 1)d 9 

G{C,i;T) = A{l+Tf £ J]^ ((1 + T) A < - l) ' 

ni,...,?i g =0 i=l 
(JV-l)d g r n, , . 

= A(i + rf ^ fj a?* 53 ( ?* ) (i + 

m,...,n g =o i=l L fc i= o ^ l ' 

(N-l)d g (jV-l)d 

^a y (i+r)^ A n e (-ir-v*^) ( m ° dj )- 

Applying Q to each side we get 

(N-l)d g 

F(C,Z;T) = A Y ( 1 + T ) M{P+k - ) tl B k i ,(N-l) d mi)) (modT N ). □ 

fcl,...,fcg=0 1=1 

Corollary 4.2. The power series F(C,c;T) has coefficients in Z p . 
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Proof. Let £ be a non-trivial element of X(c) and write C = C(/3; Ai, . . . , X g ). For i = 1, . . . , g 
we know 1 — £(Aj) is divisible only by primes above c and hence is invertible in Z p [£]. Therefore 
the polynomials Fn(C,^;T) have coefficients in Z p [£], and it follows from Galois theory (see 
also Subsection I5.3p that 

F N (C,c;T) := £ F N (C,&T) 

has coefficients in Z p . By the theorem we have F(C,c;T) = Fn{C,c;T) (mod T N ) for all 
> 0, and therefore F(C,c;T) has coefficients in Z p . □ 

since F m (a _1 ,c;T) is the sum of the F(Ci, c; T)'s where {C\, . . . , C m } is a c-admissible cone 
decomposition of o modulo m, it follows from the corollary that F m (a _1 , c; T) has coefficients in 
Z p . In particular, F m (a _1 , c; T) defines a Z p -valued measure, which we will denote /i p 'Sn- To be 
able to apply Theorem 12.91 to this measure we need to prove that Hp* m has support in 1 + qZ p . 
We will actually prove a stronger statement that will be useful later. Let Qoo be the cyclotomic 
Z p -extension of Q. Denote by Q := Q, Qi, Qa, ■ ■ ■ the subfields of Qoo/Q with Gal(Q„/Q) ~ 
Z/p n Z. Let be the cyclotomic Z p -extension of E and let E$ := E, E%, E2, ■ ■ ■ be the 
subfields of E^/E with G&\(E n /E) ~ Z/p n Z. Define integers m , mi > by Q mo = E n Qoo 
and Q mo + mi = E(m) n Qoo- By construction E mi = EQ mQ+mi is the intersection of E(m) 
and Eqq. The commutative diagram 

Cl m (^) — (1 + qp mo Z)/(l + gp m »+ mi Z) 

I (42) 
Gd(E mi /E) ► Gal(Q mo+mi /Q mo ) 

comes from Class Field Theory, where the bottom map is the restriction map, the top map 
is induced by the map h-» (Af(a)), and the vertical maps are the respective Artin maps. 
Define e > 1 to be the largest integer such that W p e c E(W g ). It is clear that we have 
e = rriQ + v p (q). We note in passing the lemma below, which will be useful later and which is 
a direct consequence of the diagram. 

Lemma 4.3. For any fractional ideal a of E coprime to p we have {N(a)) G 1 +p e Z p . □ 

We now prove our result on the support of the measure /U p 'jn. 
Proposition 4.4. The measure /j,p^ m has support in 1 +p e+mi Z p . 

Proof. Let C be a cone in a c-admissible cone decomposition of modulo m and let £ be a non- 
trivial element of X(c). Denote by /ic,£i respectively (jLq ^ with N > 1, the measure associated 
with the power series F(C,£;T), respectively Fn(C, £;T), and write C = C((3; X\, . . . , X g ). 
For k £ N 9 the algebraic integers /3 + k ■ X, are all congruent to 1 modulo q, and thus we have 
M((3 + k-X) = l (mod q) . It follows that M(/3 + k-X) = (A/(/3 + k ■ A)) € 1 + p e+mi Z p by the 
diagram above, since (3 + k • A G £? m . Thus, by Lemma [2?71 the measures have support 
in 1 +p e+mi Z p . The same is true for fic,^ since the measures /j,^ c converge (weakly) to ^c,£- 
The conclusion follows as well for /x p ' m , it being the sum of finitely many such measures. □ 

Remark 4.5. Replacing f by fp a for some a > 1 does not change the p-adic L-function, as 
we will see from the interpolation property (|4,5p it satisfies (and the unicity statement that 
follows from it; see Remark 14.1 1 j) . In particular, by taking a large enough we can force the 
measures /i p 'm to have support in 1 + p b Z p for b > 1 arbitrarily large. The proposition shows 
that (HI) is enough to imply that b is already quite large. 
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4.2. Construction of p-adic .L-functions. We are now ready to define p-adic L-functions. 
The first step is to define the p-adic equivalent of twisted partial functions. 

Proposition 4.6. For m an integer and s £Z p define Z p \^(a~ 1 1 c; s) to be the p-adic twisted 
partial zeta function given by 

Z^(a-\c;s) :=u;(M(a)r(M(a)) s J <t>- s (x)d^ m . (4.3) 

Then Z^(a~ l , c; s) is a continuous function on Z p , and 

Zp^m (a' 1 , q -k) = Z^cC 1 , c; -k) 
for all k G N such that k + m = (mod (f(q))- 

Remark 4.7. Using Proposition 14.41 we could restrict the domain of integration in (|4.3p to 
1 + p e+mi 7jp and then replace <j)- s (x) by x~ s . This would give a somewhat neater formula, 
and we will actually use this expression several times in the next subsection. However, from 
a computational point view it is better to express things as in (|4.3p . since that is how the 
computation will actually be done. 

Proof. This is a direct application of Theorem l2.9l (changing s to — s) with a n = Z m (a _1 , c; —n) 
and F = F m (cr\c;T), using Corollary 13.131 and Proposition 14.41 and the fact that for k G N 
with k + m = (mod 93(g)) we have 

u(Af(a)) m (Af(a)y k =Af(a)- k . □ 

Let x be a complex character on Cl m (E). Recall that we have embedded Q into C p and thus 
we can view x a ^ so as a p-adic character. Define the character k of Cl m (-E) by the composition 

K : Cl m (E) » Cl q (E) -^-> Gal(E(q)/E) > (Z/qZ) x T p . 

The first map is the natural surjection coming from the fact that q divides m by (HI). The 
second sends a class C to its Artin symbol o~c, where E(q) is the ray-class field of E modulo q. 
The next comes from the fact that E{q) contains the q-th root of unity, and thus we can 
associate with any a G Gal(E(q)/ E) a class a in (Z/qZ) x such that <r(£) = C a f° r au C £ W q . 
The last map sends a to w(a) G T p . For a fractional ideal a of E, relatively prime to m, it 
follows from the definition of the Artin map that k(o) = u(J\f(a)). As a consequence of the 
previous result and of Proposition 13.11 we recover the construction of p-adic L-functions. 

Theorem 4.8 (Deligne-Ribet, Barsky, Cassou-Nogues). Letm be an integer, and ifx 7^ n 1 ~ m 
assume further that c is such that x(c) 7^ K 1_m (c). Define a p-adic L-function by 

/ / \ \ -l h m (E) 

L ( As) := ( ^%Z! " l) E x(^)Z^(a-\ q s) (4.4) 

where the sum is over integral ideals Oj, relatively prime to m and c, representing all the classes 
ofCl m (E). Then L p r ^(x', s) is a continuous function on 7h v (respectively Z p \{l} ifx = K 1 ~ m ) 
and 

L p m 2(x;-k) = L m (x;-k) (4.5) 
for all k G N such that k + m = (mod f(q))- 

Remark 4.9. Assume x^™ 1 " 1 ls n °t the trivial character and let N be the Galois closure of 
E(m)/Q. Let a G Gal(E(m)/E) be such that XK m ^(a) ^ 10 Lift a to an arbitrary element 
a of Gal(iV/Q). By the theorem of Cebotarev [Ml Chap. VIII, Th. 13.4] there exists a positive 
proportion of prime ideals of N whose Frobenius in iV/Q is a. Let <t be one of these prime 

AVe see X Rm_1 as a character on Gal(_B(m)/i5) via the Artin map. 
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ideals with £ coprime to fZjv- Then c := £ nZg is a prime ideal of E, coprime to f, of residual 
degree 1, such that x K?n HO 7^ 1- 

Proof. By the previous result the sum in the RHS of (|4.4p is a continuous function on Z p . We 
now look at the factor before the sum. It is continuous at s E Z p unless x(c) = w(c) m-1 (c) a . 

Since (c) has infinite order, this can happen only for s = 1. Thus, L p ^(%;s) is continuous 
on Z p \ {1}, and also at s = 1 if %(c) / w(c) 1_m . We get the first result by noting that 
uj(c) = k(c). For the second, let k £ N be such that k + m = (mod g). The interpolation 
property follows from Proposition 14.61 Proposition 13.11 and the fact that 

m-l/„\fc-l Av / • ^ 



uj{c) m - l {c) 

Remark 4.10. Assume y^" 1-1 is non-trivial but Y^ m_1 (c) = 1- Then (|4.4p still defines a 
continuous function on Z p \ {1}. Moreover, that function is equal to L p ^(x; s) if s 7^ 1, and 
therefore, by continuity, it converges to L p ^(x; 1) as s tends to 1. 

Remark 4.11. Assume p is odd. Then the set {— k with k £ N and A; + m = (mod p — 1)} 
is dense in Z p , and there is a unique continuous p-adic function satisfying (|4.5|) . This proves 

that L p ^i(x; s) does not depend on the choice of c. For p = 2 the closure of {—k with k € N 
and + m = (mod 2)} is either 2Z2 or 1 + 2Z2. Thus (|4.5p is not enough to prove unicity. 
In that case we can use Proposition 14. 131 below!^ to conclude that L^^ix'i s ) does not depend 
on the choice of c. 

There are actually f>(q) twisted partial zeta functions or L-functions defined by these two 
results, depending on the choice of the congruence class of m modulo <p(q). However, when x 
is the trivial character we see that only for m = 1 (mod f(q)) does the corresponding p-adic 
zeta function possibly have a pole at s = 1. (See [TT] for the computation of the residue; the 
fact that it is non-zero is equivalent to the Leopoldt conjecture.) Therefore, the p-adic L- 
function, denoted simply L p m (x; s), is defined to be the function corresponding to the choice 
m = 1 (mod 99(g)) . From now on we will focus uniquely on that case, dropping the exponent 

in the notation and writing Z P)Xn (a~^, c; s) instead of Z^a^oT 1 , c; s), and so on. We will see 
below (see Proposition I4.14p that the different L-functions for various m can be recovered 
from L PiW (x; s) by twisting the character x by some appropriate power of k. 

4.3. Some properties of p-adic .L-functions. In this subsection we prove some well-known 
results about p-adic L-functions that will be useful later. The first is a direct consequence of 
Proposition 14.61 (and the remark that follows it). 

Proposition 4.12. The Mahler expansion of the p-adic twisted partial zeta function is 

Z p ^a-\v )S )=u{M{a))Y, I {{M{a))x- l -l) n d^J S ). □ 

n>0 Jl +P e+miZ P 

From this, we deduce the analyticity of p-adic L-functions. 

Proposition 4.13. Let B e be the open ball in C p of center and radius p e ~ l /^ p ~ l \ Then the 
p-adic L-function L Pjm (x; s) can be extended to an analytic function on B e , if x is non-trivial, 
and to a meromorphic function on B e with a pole of order at most 1 at s = 1, if x is trivial. 

Proof. We first prove that the p-adic twisted partial zeta functions Z Pim (o _1 , c; s) can be 
extended to analytic functions of radius^ p e-1// ( p_1 ). For x £ 1 +p e+m-1 Z p it follows from 



^The proposition only applies to m = 1 (mod 2), but it is straightforward to generalize. 

^We say that an analytic function has radius r if it converges on the open ball in C p of center and radius r. 
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Lemma 14.31 that |((AA(a)x x ) — l)™| < p en , and since the measure [i p ' m is of norm < 1 

and -F m (a _1 , c; T) has coefficients in Z, p we conclude by the previous proposition that the n-th 
Malher coefficient of 2p ]m (a _1 , c; s) has p-adic absolute value < p~ en . The result then follows 
from Corollaire 2(c) of [I, p. 162]. (See also Theorem 15.91 ) 

Denoting by g(s) the inverse of the factor before the sum in the RHS of (|4.4p we have 

g(s) := x(c)(c) 1 ~ s - 1 = x(c) exp p ((l - a) log p (c» - 1 



where exp p and log p are respectively the p-adic exponential and logarithm functions (see [2 
§5.4]). The function g(s) is analytic on B e since | log p (c)| < q~ e , using again Lemma fOl and 
the fact that the p-adic exponential function has radius p -1 /^- 1 ). Furthermore, from the 
properties of the p-adic exponential and logarithm functions we see that g has a simple zero 
at s = 1 if x( c ) = 1) an d does not vanish otherwise. Now if x is non-trivial we can proceed 
as in Remark 14.91 and choose c such that x( c ) 7^ lj this proves the result for the first case. 
Otherwise we write g{s) = (s — l)h(s) where h(s) is an analytic function non-vanishing on B e 
and the second case follows. □ 

The next result establishes that the p-adic L-function L p "^(x;s) for any m in Z can be 
recovered from the p-adic L-function (corresponding to m = 1). 

Proposition 4.14. For any integer m and any s G Z p; assuming s / 1 if \ = K l ~ m , we have 

L { p m 2( X ;s) = L p , m ( X K 1 - m ;s). 

Proof. Let a be an integral ideal coprime to p. Then Z p r ^!(a~ 1 , c; s) = K(a) m-1 i? Pjm (a -1 , c; s), 
since k(o) = uj(M(o)), and substitution in (|4.4p yields the result. □ 



A key property of p-adic L- functions is that they are Iwasawa analytic functions (see [2"5]): 
we will show this in the proof of the next theorem. Let E x be the subextension of E(m)/E 
fixed by the kernel of X - After Greenberg [18] we say that x is °f type W if E x C K^J 10 I 



Theorem 4.15 (Deligne-Ribet). Fix a topological generator u of 1 + p e 7L p . Then there exists 
a unique power series 3p, m (x; X) in Q p (x)[[X]\ — or X^ 1 Q p (x)[[X]] if x is trivial — called the 
Iwasawa power series of x {with respect to u) such that, for all s G Z p (excepting s = if x is 
trivial) we have 

L p , m (x;l-s) = 3 p ^{ X ;u s -I). (4.6) 

Moreover, if x is trivial or not of type W then 3p, m (x; X) has coefficients in Z p [x]. Otherwise 
there exists a non-trivial root of unity £ of order dividing p mi such that (£(1 + X) — l)3 PtVa (x] X) 
has coefficients in 7L p \x\. 

Proof. Unicity is clear by fj4.6|) since the set {u s — 1 with s G Z p } = p e 7L p admits as a limit 
point. In particular, if the Iwasawa power series exits it does not depend on the choice of c. 
We will use this fact below by choosing prime ideals c satisfying additional properties. We 
now prove existence. For x G Z p with (x) G 1 + p e 7L p define 

log„n 



so that 



^^See the discussion after Corollary 14.21 for the notation used in the theorem and its proof. 
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In particular, we can use this equation with x := N(a), for some ideal a coprime to p, by 
Lemma [4T3l We define three power series, with coefficients in Z p , Z p , and Z p [x] respectively, 
as follows. 

N(a; A) := M{a) £ f~ C ^(a))\ %l = ^ (fl)(1 + X) -UWW) , 
,4(a,c;A) := W x-'^fAd^X^ ! x-\l + Xf^d^ m , 

C(c, X ; X) := X (c) £ ( £ ^ (c) ) X* - 1 = X (c)(l + X)^ - 1. 

Then for all s G Z p , using (j4.3|) for the first equality, we have 

iV(a; u s - l)A(a, c; u s - 1) = ^, m (o _1 , c; 1 - s), 

C7(c, X ;^-i) = x(c)(c) s -i. 

We define 

h m (E) 

3 p , m (x;X):=C(c, X ;X)- 1 £ x(^ 1 )N(a l ;X)A(a l ,c;X) (4.7) 

i=i 

where the sum is over integral ideals Oj, relatively prime to c, representing all the classes of 
Cl m (E). For all s G Z p such that C(c, x;^ s — 1) ^ it follows from (|4.4p and the equalities 
above that (|4.6f) holds. 

We now consider several cases, not necessarily disjoint. Assume first that the order of x 1S 
not a power of p. We reason as in Remark 14.91 with a € Gal(E(m)/E) such that the order of 
x(&) is not a power of p. Then x(c) — 1, the constant coefficient of C(c, X) -X"), is a p-adic unit. 
Therefore C(c, x; A) does not vanish on p e Z p and is invertible in Z p [x][[X]]. This proves that 
^p,m(x;-X) ^ Z p [x][[-X]] and satisfies (|4.6|) for all s G Z p . 

Assume now that x is such that £7 X n E^ = E. Let N be the Galois closure over Q of 
the compositum of E x and E\. Then there exists a G Gal(A/Q) such that <?\e x is trivial 
but is non-trivial. Let £ be a prime ideal of N, coprime to fZjv, whose Frobenius is 
equal to a, and let c := £ n Zg. Then c has residual degree 1, and x(c) = 1 so the constant 
term of C(c, Xi-%) is zero. Also, by construction we have C u (c) G Z p , and since this is the 
coefficient of X in C(c,x;X) it follows that C(c,x; A) = A*7(X), where U{X) G Z p [x][[A]] 
is an invertible power series. Therefore C(c, x; A) vanishes only at 0, that is, for s = 0, and 
3 p , m (x; A) G A _1 Z p [x][[A]]. This proves the result when x is trivial. When x is non-trivial, 
the limit of 3 p>m (x; A) when A — >■ exists and is finite^ Therefore the coefficient of A -1 in 
^p,m(x; -X) is zero. This proves the result in the non-trivial case. 

Before looking at the other cases we need additional notation and results. Recall that mo 
and mi are such that Q mo = E f] Qoo and Q mo +mi = E(m) n Qoq. Let £o be a root of unity 
of order p mi . Define a function on 1 + qp m °Z p by 

This map is a group homomorphism with kernel l + gp mo+mi Z p . Composing with the function 
on the top in diagram (j4.2j) we get a character 

Co 

on the ray-class group Cl m (E). It is easy to see from its construction that this character 
generates the subgroup of characters of Cl m (E) of type W. Since e = uiq + v p (q), we have 



ii 



It is equal to L p , m (x; !)• 
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log u/(qp m O) 

logp u/ (qp m °) G Z* and without loss of generality we can replace £o by £ p to get the 

character p : a i— > ^jf 1 *^ ^, which still generates the group of characters of Cl m (i?) of type W. 
For wGZwe compute 

A(a,c;^(l + X) - 1) = / x-\Q{l + X)f^ dp, a p % = A(a,c;X). 

We also have N(a;Q(l + X) - 1) = p~ v (a)N(a;X) and, for any character if) of Cl m (E), 
C(c,V;$(l - 1) = C(c,ipp v ;X). We conclude that 

3 p , m (if) P v ;X) = 3 p , m (il>;®{l + X) - 1). (4.8) 

Assume now that \ is OI tyP e W. Then x = P v f°r some v G Z. Using (j4.8|) with ^ the 
trivial character, we find that (£q(1 + X) — l)3 Pjtn (x; .X") ^ ^p[x][[-^]]> which proves the result 
in this case, taking £ := £q. 

Finally we consider the case E x <f_ E^, in which \ is not of type W and, in particular, 
E x D -Eoo ^ E. In. this case we can write x = ^p 11 for some non-trivial character i/> of Cl m (E) 
satisfying E^, n -Eoo = -E 1 and some v G Z. Since the Iwasawa power series for ?/; is in Z p [x] [[X]] 
by the argument above, it follows from (|4.8p that the same is true for 3 p ^ m {x',X). □ 

As a first application we use this result to bound the size of the values of p-adic L- functions. 

Corollary 4.16. Let s G Z p , with s^l i/% w trivial. Then 

1 if X is n °t °f tyP e W > 

p-V(P _1 ) ij^- i s of type W and non-trivial, 



\L P} m(x;s)\ p < < 



p e | 1 — sL 1 if x is trivial. 



ip 

Proof. The result is clear if x is not of type W0 Assume x is of type W and non-trivial. In 
the notation of the theorem we have u p (£ — 1) < l/(p — 1) and v p {u l ~ s — 1) > v p (q). From 
the identity ab — 1 = (a — 1)6 + 6 — 1 we get f p (£tt 1_s — 1) < l/(p — 1) and the result follows. 
For trivial x we have i> p (u 1_s — 1) = u p (l — s) + e, which proves the result. □ 

5. Computational methods 

In this last section we show how to use the results of the previous sections to compute values 
and representations of p-adic L-functions explicitly. Note that all computations will involve 
only p-adic integers — approximated by integers as we explain below — and that we will need 
to deal with p-adic rational numbers only in the last subsection (and we will do it somewhat 
indirectly) . 

In discussing the computation of p-adic approximations we will follow certain terminological 
conventions. Let M > 1 be an integer. For a G Z p we denote by a mod p M the unique integer 
a with < a < p M such that \a — a\ p < p~ M . By "computing a to the precision p M " 
we mean computing a mod p M . Let L be a finite dimensional Z p -lattice with v\, . . . , v n a 
(fixed) basis of L. For a G L, by "computing a to the precision p M (with respect to the 
basis v±, . . . , v n )" we mean computing the p-adic numbers a%, . . . , a n to the precision p M , 
where a = a%v% + • • • + a n v n . Note that in what follows the basis v\, . . . , v n is usually not 
stated explicitly but it should be clear from the context what it is. Let N > 1 be an integer 
and let F(T) G Z P [[T]] be a power series. By "computing F to the precision (p M ,T N )" we 
mean computing the first N coefficients of F to the precision p M . Let p G .A4(Z p ,Z p ) be a 
measure with values in Z p . Its associated power series Fu(T) must therefore lie in Z p [[T]]. 
By "computing p to the precision (p M ,T N )" we mean computing the power series F^ to the 
precision (p ,T N ). Finally, let / be a continuous function in C(Z p ,Z p ). Then its Mahler 



12 In this case we also get the well-known statement: L p>m {x\ s ) = £p,m(xi 0) (mod p e ). 
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coefficients (f n )n>o are all p-adic integers and tend p-adically to zero. For M > 1 we denote 
by Nf(M) the smallest integer N > such that \f n \p 5; P~ M for all n> N. By "computing / 
to the precision p M " we mean finding an integer N > Nf(M) and computing the coefficients 
/o) - j /jV-l t° the precision p M . 

In the complexity estimates below it is assumed that fast multiplication algorithms are 
used. Therefore, for example, it takes 0~(M log p) bit operations to multiply two rational 
p-adic integers to the precision p M , and it takes 0~(NM logp) bit operations to multiply two 
power series in Z P [[T]] to the precision (p M , T ). Here, to simplify the complexity expressions, 
we have used 0~ -notation: g G 0~(f) if there exists c > such that g G 0(/(log/) c ). 

Finally, we will assume that the necessary data to work in the field E have been computed. 
In particular, we assume that an integral basis, say (0\, . . . , 0^), is known. We will express the 
elements of E with respect to this basis. Also, we assume that the class group, the group of 
units, and the ray-class group modulo m are known. Algorithms to perform these tasks can 
be found in [8] and |9j; see also [3]. 

5.1. Computations with continuous functions. Let / G C(Z p ,Z p ). For 7V > 1, we com- 
pute the first N Mahler coefficients of / with the following algorithm. 

Algorithm 5.1 (Computation of Mahler coefficients). 
Input: / G C(Zp,Zp). 

Output: The first N Mahler coefficients of / to the precision p M . 

1. For n = to AT — 1, do f n -<— f(n) mod p M . 

2. For j = 1 to N — 1, for n = N — 1 to j (decreasing), do f n ^—fn — fn-i mod p M . 

3. Return / , . . . , Jn-X- 

Lemma 5.2. Assume for x G Z p that it takes 0(C) bit operations to compute f(x) to the 
precision p M . Then Alaorithm \5.1\ computes the first N Mahler coefficients of f to the precision 
p M in 0(NC + N 2 M log p) bit operations. In particular, for s G Z p it takes 0~(NM 2 log 2 p + 
N 2 Mlogp) bit operations to compute the first N Mahler coefficients of (j) s to the precision p M . 

~(i) 

Proof. Let /A denote the value of f n after j iterations of the main loop in Step 2. We claim 
for < j < N - 1 that 

f(V"/)(0) mod p M for < n < j, 

y (V j f)(n - j) mod p M for j < n < N - 1, 

where V is the finite-difference operator defined by (V/)(s) := f(s + 1) — f(s). The claim 
follows for j = by the initialization in Step 1 since V° is the identity. Assume now that the 
claim holds for some j. If < n < j then /A — fn and the result is proved. If n > j + 1 
then 

fn J+l) = fi j) ~ Bi mod p M = (VV)(n - j) - (VV)(n - j - 1) mod p M 

= (V(VV))(n -j -I) mod p M = (V j+1 f)(n - {j + 1)) mod p M 

and the result follows by induction. In particular, at the end of the algorithm we have 
fn = (V n /)(0) mod p M = f n mod p M , where (f n )n>o are the Mahler coefficients of / (see 
[26} §2.4]). This proves that the algorithm returns the correct result. We now estimate its 
complexity. The initial step takes 0(NC) bit operations by definition and the second step 
takes 0(A 2 Mlogp)El This proves the first complexity statement. 

Now we turn to the computation of 4>s(x) mod p M . We assume s is given by its approxi- 
mation s mod p M , which we will still denote s by abuse. We can also replace x by x mod p M 



fU) 

J n 



l^The cost of computing the remainder modulo p M is also 0(M log p), since — p M < f n — f n 
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without loss of generality. If p divides x then (j) s {x) = 0. We now suppose that x 6 Z*. For 
p odd, assume we have computed and stored the values uj{o) mod p M , for o = 1, . . . , p — 10 
Then we can compute oj(x) mod p A:f in 0~(M log p) bit operations, since w(x) = a; (a), where 
a := xmod p. The computation of oj(x) for p = 2 is trivial. Then (x) = x/ui(x) is computed 
to the precision p M in 0~(M log p) bit operations. Assuming fast exponentiation, it takes 
0~(M 2 log 2 p) bit operations to compute (x) s Hence it takes 0~(M 2 log 2 p) bit operations 
to compute the value of 4> s (x) to the precision p M . Combining this with the first complexity 
result completes the proof of the last statement. □ 

To use the Mahler expansion to compute values of a continuous function we need to compute 
binomials coefficients (°) to the precision p M for n = 0, . . . , N — 1. (We will also need these 
coefficients for the computation of measures and Iwasawa power series.) We use the following 
algorithm. 

Algorithm 5.3 (Computation of binomial coefficients). 
Input: sSZ,,. 

Output: The binomial coefficients for n = 0, . . . , N — 1, to the precision p M . 

0. For n = 1 to N, do v n <— v p (n) and u n <— {np~ Vn )~ 1 mod p M . 

1. Let V be the largest integer v > such that p v < N — 1. 
Set s •<— s mod p M+ ^. 

2. 8etA<-l,B<- 0, 6 <- 1. 

3. For n = 1 to iV - 1, do 

If s — n + 1 = 0, set bk 
b <— v p (s — n + 1), a <- 
A <— au n A mod p M , i? 
6 n ^- vlp B mod p M . 

4. Return &o, ■ • • , &JV-1- 

Remark 5.4. The precomputations in Step need to be done only once for fixed and M. 

to the precision p M in 0~(NM log p) 



— for k = n to N — 1 and go to Step 4. 
(s — n + l)p -6 mod p M , 

- B + b — v n , 



Lemma 5.5. Ala orithm \ 5. 3\ computes (g), Q), 
6zf operations. 



Proof. For n > and x G N define 



P 



if C) = 0, 
otherwise. 



From the recurrence relation satisfied by binomial coefficients it follows that 

xl (x — n + \)p- v p( x ~ n + l ) r x 
n np~ v p( n } n — 1 

From this one can see by induction that at the end of n-th iteration of the loop in Step 3 
we will have A = [*] mod p M and B = ))■ Thus the algorithm returns ( s ) mod p M for 
n = 0, A r — 1. Now for n > 1 we can write (*) = 5(fZi) an d therefore 



< 



n 



11 



These can be easily computed using Hensel's Lemma. 
-^For p odd, one could instead compute x t directly, with t as in the proof of Proposition 15.101 
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It follows that (1 + T)P M+V = 1 (mod p M ,T N ) and thus 

{1 + T) S = (1 + T) 5 (mod p M ,T N ). 

Therefore f * ) = f*) (mod p M ) for n = 0, . . . , N — 1, and hence the algorithm returns the 
correct result. 

Lastly we estimate the complexity of the algorithm. For an integer n > 1 we can compute 
v p (n) and np~ Vp ^ using v p (x) + 1 divisions by p. Therefore Step performs 

N + [N/p\ + [N/p 2 \ + • ■ • < N/(p - 1) 

divisions and therefore takes 0~(NM logp) bit operations. However, if we use the same method 
in Step 3 we may end up needing many divisions if s — n + 1 has a very large p-adic valuation. 
A better way to proceed is to use a divide-and-conquer algorithm, so that the computation 
of a and b can be done in 0(log(M + V)) divisions. The computation of A and b n takes 
(X(Mlogp) bit operations, the computation of B is negligible, and, since V G 0(log(A)), 
Step 3 takes 0~(NM logp) bit operations. This completes the proof. □ 

Remark 5.6. The algorithm can be improved in the following way. In Step and Step 3 
we can keep a counter, say ct, that we set to p the first time we encounter an integer with 
a non-zero p-adic valuation. Then at each step we decrease ct by 1. If the value of ct is 
non-zero then the p-adic valuation of the corresponding integer is zero. Otherwise we compute 
the valuation using the method explained above and reset ct to the value p. This gains a 
factor p in computing the p-adic valuation and the prime-to-p part. It does not change the 
total computation cost estimate however. 

Once we have computed sufficiently many Mahler coefficients of / to the precision p M we 
can use the algorithm below to compute values of /. 

Lemma 5.7. Let f G C(Z p ,Z p ). Assume we have computed f to the precision p M with respect 
to N > Nf(M). Then for all s G Z p we can compute f(s) to the precision p M in 0~(NM logp) 
bit operations. 

Proof. Let /o, . . . , fN-i be the first ./V Mahler coefficients of / to the precision p M . Then 

N-l , s 

n=0 

The binomial coefficients are computed using Algorithm 15.31 in 0~{N M\ogp) bit operations, 
and the computation of the sum also takes 0~(./V"M logp) bit operations. □ 

Another reason to compute Mahler coefficients is for approximating integrals. 

Lemma 5.8. Let f G C(Z p ,Z p ) and let [i G A / f(Z p ,Z p ). Assume we have computed f to the 
precision p M with respect to N > NAM) and have computed fi to the precision (p M ,T N ). 
Then we can compute f fdfi to the precision p M in 0~(NM\ogp) bit operations. 

Proof. Write Fq, . . . , Fn-i (respectively /o, . . . , /jv-i) f° r the first coefficients of F^(T) 
(respectively Mahler coefficients of /) computed to the precision p M . We have 

N-l 

n=0 

This computation takes 0~(AMlogp) bit operations. □ 



ffdfx=J2fnFn (mod^ 

" : n 



From these results it is obvious that having the best possible upper bounds on Nf is crucial 
for getting the best complexity estimates. In the next subsection we consider this problem for 
the functions that interest us. 
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5.2. Analyticity and Mahler coefficients. A power series in C P [[X]] is restricted if its 
coefficients tend to zero or, equivalently, if it converges on O p := {x 6 C p such that \x\ p < 1}. 
Let / : r L p — > C p be a function. We say / is analytic if there exists a restricted power series 
f(X) € C p [[X]] such that 

f(x) = f(x) for all x 6 Z p . 

For h> we say / is locally analytic of order h if there exist restricted power series f a ,h{X), 
with < a < p — 1, such that 

/(x) = /o,fc((a: - for all x £ a + p h Z p . (5.1) 

Note that the series f a h are uniquely defined. An analytic function is therefore a locally 
analytic function of order and one can verify that, if / is analytic of order ho, then it is 
analytic of order h for any h > Hq. Clearly a locally analytic function is continuous. We 
will see below that the fact that a continuous function is locally analytic has some important 
consequences for the rate of convergence to zero of its Mahler coefficients The norm of a 
restricted power series /, denoted ||/||p jQ o, is defined as the maximum of the absolute values 
of its coefficients. It can be computed thanks to the following (see \26\ Prop. 1, §6.1.4]): 

11/llp.oo = max \f(x)\ p . (5.2) 

Let / be a locally analytic function of order h. We define the /i-norm of / by 

M h (f) := max ||/ a ,h[Loo- 

0<a<p h -l 

It follows from <&T§ and ([O]) that 

M h (f)= max \f(a + p h x)\ p . (5.3) 

0<a<p h -l 

xeOp 

Theorem 5.9 (Amice). Let f be a function in C(Z p ,C p ) with Mahler coefficients (f n )n>o- 
Then f is locally analytic of order h if and only if 



f n 



-»• 0. 

v 



\njp h \\ 

Moreover, if f is locally analytic of order h then for all n > we have 

\f n \ P <M h (f)\[n/p h \\\ p . 

Proof. The first assertion is Corollaire 2(b) of [U p. 162]. For the second assertion we use 
the notation and results of Corollaire 1 of [TJ p. 157]. Since the sequence u n = n is very well 
distributed (tres bien repartie), it follows from Corollaire 1(c) that 

h 

v p (n\/s n , h ) = J>//J " ELV/J = v P ([n/p h \l). 

k>l k=l 

Using Corollaire 1(b), we find that (L n /^' l J'(n))n>o ^ s a norma l basis (base normale) of the 
Banach space of locally analytic functions of order h. Hence Mh(f) = sup n>0 \ fn/\n/p h \ \\ p , 
and the result follows. □ 



l^The main reference for these results is the article of Amice [I]; see also [10| for a more accessible 
presentation. 
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We are interested in finding optimal upper bounds for N ( p a . We remark that the function 
j is locally analytic of order 1 if p is odd and of order 2 if p = 2. Indeed, for x G a + qZ p 



with a G Z* we have 

W . = (a) .(£z£ + 1 )' = <a> .g(^(£z£)". (5 ,, 

In this way Theorem 15.91 provides bounds on the Mahler coefficients of 4> s . But we can do 
better with the following result. 

Proposition 5.10. Fix s G Z p and let (4> n )n>o be the Mahler coefficients of 4> s . Then for all 
n > we have 



Kip < 



VK2J+1 if p = 2, 
\n\\ p if p is odd. 



Remark 5.11. This result is close to optimal for odd p. Indeed, it implies that \<p n /n\\ p < 1 
for all n > 0. On the other hand, we know this quantity is also bounded from below; otherwise 
Theorem 15.91 would imply that (f) s is analytic, and in general it is not. 

Proof. Assume p is odd and let B be a positive integer. By the Chinese Remainder Theorem 
we can find a positive integer t such that 

t = s (mod p B ), 
t = (mod p — 1), 
t > B. 

Then x l = uo(x) t (x) t = (x) s (mod p B ) for x G and x l = (mod p B ) for x G pL v . Hence 
~ x^\p < p B for all x G Z p . It follows from (|2.4p and Theorem [5^9] that 

\tp n \ p < max{|n!| p , p~ B ] 

and we obtain the result by taking B large enough. 

For the case p = 2 a similar proof works, provided s is even. But for odd s we need to use 
another approach. So we assume s G 1 + 2Z2. Let B and t be positive integers such that s = t 
(mod 2 B ). As above we have 

\(f )s {x)-uj{x)x t \ 2 <2- B (5.5) 

where uj{x) := if x G 2Z2. We turn now to the computation of bounds on the Mahler 
coefficients (a n ) n >o of w(x)x . Let i be a fixed square root of —1 in Q2. Then x \-t (±i) x are 
continuous function^ on Z2 and 

= ^(Hr -n = \ I> - i) n (^ n - 1) (*) • 

n>0 ^ ' 

Thus the Mahler coefficients (w n ) n >o of w satisfy 

+00 if n = (mod 4), 

V2(wn) = { n/2 if re = 2 (mod 4), 

k (n-l)/2 if re ee 1,3 (mod 4). 

In particular, V2{w n ) > [ n /2j for all n > 0. We now apply the following lemma. 



But they are not analytic functions. 
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Lemma 5.12. Let f and g be two continuous functions with Mahler coefficients (f n )n>o an d 
(<7n)n>o- Then the Mahler coefficients (c„) n >o of fg satisfy 

v p {c n ) > min (vp(fk) + min v p (g n 
0<k<n \ n—k<m<n 

Proof of the lemma. Let k, m > be two integers. Then 

k+m 



n=max(fc,m) 



for some u n (k,iri) G Z. Hence 
f{x)g{x) 



E ^fc) (™) = E (l> E ^(*.™))(*) 

and the result follows from the expression of c n that can be derived from this equality. □ 
We take f(x) = oj{x) and g(x) = x l . Then 

V2{a n )> min (i>2(wfc) + min v 2 (m!)) = min (^(^A;) + V2((n — k)\) ) 
> min (Lfc/2J + L(n - fc)/2j) > [ri/2j - 1 

0<fc<n 

and the result follows from this estimate, taking B sufficiently large in (|5.5p as before. □ 
Corollary 5.13. For every positive integer M we have 



2M + 2 ifp = 2, 
pM if p is odd. 



Proof. The result is clear for p = 2. For p odd, it is enough to prove that the p-adic valuation 

of (pM)\ is at least M. But v p ((pM)\) = £ [pM/p fe J > M, and the result follows. □ 

fc>i 

Recall that for x S Z* with (x) G 1 + p e Z p we have set 



£ u (x) :-- 



logpM 

where it is a fixed topological generator of l+p e 7j p . For an integer I > we define a continuous 
function tp£ in C(Z p ,Z p ) by 



( x -i ir ( .r}ei + ,/^. 



otherwise. 



We have aT^l + SJ^O") = E^>o^(^)^ if <») G 1 + P e % P - These functions appear in the 
construction of the Iwasawa power series and will play an important part in their computations. 

Proposition 5.14. The Mahler coefficients (tp£ n )n>o of ip£ satisfy 
for all n > 0. 

Proof. It is clear that the function ip£ is locally analytic of order e. The e-norm is l/|£!| p , by 
(15.31). The result follows from Theorem 15^91 □ 
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Corollary 5.15. For every positive integer M we have 

N^(M)<p e (pM + £). 

Proof. It is enough to prove that v p ((pM + ■£)!) > M + v p (£l). But this is clear from the 
facts that v p ((a + 6)!) > v p (al) + v p (bl) for any two non-negative integers a and b, and that 
v p ((pM)\) > M. □ 

Remark 5.16. This upper bound is in general quite far from optimal. For example, with 
p = 3, e = 1, and £ = 10 the corollary gives an upper bound of 210 for M = 20, but 
computations give A^ 10 (20) = 85. As the complexity of the computation of Iwasawa power 
series depends heavily upon the size of N^ t (M) — see for example Theorem 15.241 — it is a 
good idea to precompute the values N^ e (M) for small i and M and to use these instead in 
the computations. 

5.3. Computation of p-adic cone zeta functions. From (|3.16|) we see that, once a cone 
decomposition has been computed^ the computation of the p-adic twisted partial zeta func- 
tions, and in turn that of the p-adic L-functions, boils down to the computation of p-adic 
cone zeta functions. We now turn to this computation, but first explaining how to deal with 
the several zeta functions associated with different additive characters all at once. Define the 
etale algebra 

TZ:=Q P [X]/(X C - 1 + + 

and an additive character 5 from %e to 1Z by setting 

E(a) := rf (5.6) 

for each a £ Ze, where rj is the image of X in 7Z and a is any positive integer such that a = a 
(mod c). The next result is straightforward. 

Lemma 5.17. Let T-r be the trace oflZ/Q p and let a £ Z#. Then 

T 1l (E(a))= ^ £(a). □ 

Note that is trivial to compute since it is Q p -linear and we have 

m , r,-, I c — 1 if c I a, 
T n (r, a ) = { 1 ' 

I —1 otherwise. 

Let C := C(/3; Ai, . . . , \ g ) be a c-admissible cone. For JV > 1 we define 



F N (C,c;T) :=T n 



(JV-l)d 



A(C,~) £ (l + T)^-A)n^,(A f _i),(S(A J )) 



i=i 



(5.7) 



where A(C, H) := S( ( 9)/rj| =1 (l — H(Aj)) and is extended in the natural way to 7£[[T]]. 
It follows from Lemma [5T7] and Theorem O that Fat(C, c;T) = F(C,c;T) (mod T^). We 
will use the expression above to compute approximations of F(C, c; T) and of its associated 
measure n c p c . We define 

Z P (C,c;s) := J ^ s (x)d^ c . (5.8) 

We now determine some computation costs. These results, or at least their proofs, will be 
useful later to estimate the complexity of the computation of p-adic L-functions; see Subsec- 
tion [53J The first step is the computation of values of the rational functions B^r. 



l^This will be the topic of the next subsection. 
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Proposition 5.18. If K is a non-negative integer then 

K 



Bq,k{x) 



X 



x-1 

and for < k < K we have the recurrence formula 



x + l 



B, 



k+l,K 



» 



-iy 



K 



+ B k . 



K 



k + l J \x-l 

Proof. We first establish another expression for the B^ #-(a;)'s. 

Lemma 5.19. For k > let Coefffc denote the linear map that sends a polynomial in Q(x)[X] 
to the coefficient of its monomial of degree k. Then 



Bk,K(x) 



Proof of the lemma. We compute 

N k k . 
-x \ / n 

^1 ) ^ \k 

7 n=k v 



-X 



Bk,x(x) 



X 



CoefTi 



n—k 



+ x 



K+l 



1 



+ x 



1 



and the conclusion follows by evaluating the sum. 
Now define 



A K {X) :-- 



1 



+ X 



K+l 



k K 

Coeff A 

n=k 



+ x 



+ x 



□ 




and let a k denote the coefficient of X k in Ak{X). Then Bq^k = Ak{0), which gives the first 
assertion. Since 

1 „\ lx 



A K (X) 



we see that 



1 



+ a k 



+ X 



1 



K+l 

+ x) - J 



K + l 

k + l 



K-k 



□ 



for < k < K. The second assertion follows by induction. 

Corollary 5.20. Let a G 7L>e, coprime with c, and let K > 1 be an integer. Then we can 
compute Bq^IR (a)), . . . , Bx^(S(a)) G TZ to the precision p M in 0~(KMclogp + dlogc) 
binary operations. 



,K+1\ 
\K+l) 
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Proof. We apply Lemma 15.51 to precompute the binomial coefficients ( J" ) , . . . , 
0~(KM log p) bit operations. Since the prime ideal c is of degree 1 there exist integers 
ti, . . . , td G {0, . . . , c — 1} such that 0j = t{ (mod c) for i = 1, . . . , d. We assume ti, . . . , td 
have been precomputed (using, say, [HI Algo. 1.4.12]). Write a := a\6\ + ■■■ + a^Qd G Zg. 
Then a = a (mod c), where a := a±t± + ■■■ + adtd mod c is computed in CT(dlogc) bit 
operations. Next we precompute {n a /{n a — 1)) K in 0~(Mc log K log p) bit operations. From 
this we get i?o,i^( r / a ) at negligible additional cost. Then, using the induction formula and the 
precomputed values, the computation of -B/ c+ i j x(?7 a ) from Bf t ^{ r l a ) takes only 0~(Mclogp) 
bit operations. □ 

We now can give our first estimate. As we want our results to be valid for several different 
cones at once, we will express these estimates using d and not g (the number of generators), 
using the fact that g < d. 
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Theorem 5.21. For any positive integers M and N we can compute the measure fi pC to the 
precision (p M ,T N ) in 0~(N d+1 d d+3 Mclogp) bit operations. 

Proof. We compute F(C;T) using (|5.7|) . Applying the previous proposition we precompute 
B k ( N _ 1 \ d (S(X i )), for k = 0, . . . , (N — l)d and i = 1, . . . , d, in 0~{Nd 2 Mc\ogp) binary 
operations. With these values precomputed each computation of the inner product in (|5.7p 
takes 0~{dMc\ogp) bit operations. Now let a := M(f3 + k ■ A). We compute (1 + T) a to the 
precision (p M ,T N ), using Lemma l5.5l in 0~{N MXogp) bit operations, after having computed 
a to the precision of p M+v (in the notation of Lemma f5.5p . The main part of the computation 
of a is the computation of the norm, which boils down to the computation modulo p M+v of the 
determinant of a dx d matrix; see [U §4.3]. This takes 0~(d s (M + V) \ogp) bit operations. The 
multiplication of the power series with the inner product takes O'(NMclogp) bit operations. 
The sum has 0(N d d d ) terms so the computation of the sum requires 0~(N d d d (N + d)Mc\ogp) 
bit operations. The multiplication by A(C, 5) and the computation of the trace take negligible 
time compared to the computation of the sum. The conclusion follows by putting everything 
together and simplifying. □ 

Remark 5.22. It is a good idea to retain the values of Bj. t K(ri a ) in order to reuse them if 
several generators have the same image under H. This does not affect the estimate of the cost 
of the computation since that estimate is dominated by the cost of computing the sum. 

Once the measure fj, Pi c has been computed we can use it to compute values of Z p {C, c; s). 

Corollary 5.23. With a precomputation of cost 0~(p d+l d d+s M d+2 c) bit operations, for any 
given s in Z p we can compute Z p (C, c; s) to the precision p M in 0~(p 2 M 3 ) bit operations. 

Proof. We use the integral expression (|5.8p for Z p (C, c; s). For this we need to compute /j, c c 
to the precision (p M ,T N ) with N > N f (<f> s ). By Corollary ETJ we can take N = pM + 2; 
the cost of the precomputation comes from the previous theorem. To perform the integration 
we need to compute the first N (= pM) Mahler coefficients of <j) s \ by Lemma 15.21 the cost is 
0~(p 2 M 3 ). The cost of computing the integral, as given by Lemma 15.81 is 0~{pA / I 2 ). □ 

Recall that e is the largest positive integer such that W p e C E(W q ) and that u is a fixed 
topological generator of l+p e Z p . Denote by J P (C, c; X) the Iwasawa power series of Z P (C, c; s), 
that is, the (unique) power series in Z p [[X]] such that Z p (C, c; 1 — s) = J P (C, c; u s — 1) for all 
s G Z p . It is easy to adapt the proof of Theorem 14.151 to show that this power series exists 
We now give an explicit formula for 3 P (C, c; X) modulo (p M ,X L ). 

Theorem 5.24. Let L and M be positive integers and let K = (p e (pM + L) — l)d. Then 
3 P (C, c; X) 

■ K g 



(mod p M ,X L ^ 



A(C, E) *f(P + k- A) _1 (l + Xf^+^ [] B ki>K (E(Aj)) 

fci,...,fc 9 =0 i=l 

Hence we can compute 3 P (C, c; X) to the precision (p M ,X L ) in 0~(p ed d d+3 (pM + L) d M 2 Lc) 
bit operations. 

Proof. We start with a useful lemma. 

Lemma 5.25. Let f be a continuous function on 7L p with Mahler coefficients (f n ) n >o- Let M 
and N be integers with M > 1 and N > Nf(M). Let fx be a measure of norm < 1. Assume 
there exist a finite set A of elements of Z p and an element c a in O p for each a £ A such that 

F^{T) = Y J ^ + T) a (modT N ). 



-^Actually, the definition of 3 P (C, c; X) is given in the proof of the next theorem. 
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Then 



fdfi 



£ 

a£A 



Caf(a) 



< p 



—M 



Proof of the lemma. For each a£Z p let <5 a be the Dirac measure at a, and define the measure 
p:= (j,- J2aeA c ^cl- Then 

9 dp = / gdfj,-^ c a9(a) 
J aeA 

for all (7 G C(Z p ,C p ). From this and the facts that < 1 and |c a | p 



that /i has norm < 1. By Lemma 12, 71 the associated power series is divisible by T , and thus 



< 1 for all a £ A we see 

N 



f dp 



n>7V 



< sup |/ n L < p 



-M 



□ 



n>N 



We apply the lemma repeatedly with \i = fi c c and f = ipi for each £ = 0, . . . , L and using 



(|5.7p for the definitions of the set „4 and the coefficients c a 
N = p e (pM + L). We have 



By Corollary 15,151 we can take 



J p (C,c;X) = [ Mx)d^ c X l = E / Mx)d^,cX £ (mod X 



= T n 



L-l 



(JV-l)d 



& ■ a)) n Ski, 



,(N-l)d 



L-l 



(mod p M ,X L ) 



the expression for 3 P (C, c;X) coming from the fact that ^ ipi(x)X = 
(mod X L ). ' £=0 

We now estimate the cost of computing 3 P (C, c; X) by this formula. As above, precomputa- 
tion of the B^^s has a cost of 0~(p e+1 dM 2 Lc) binary operations, after which we can compute 
each inner product in 0~(dMc log p) bit operations. The computation of a := C u (M(f3 + k ■ A)) 
to the precision p takes 0~(M(M+e) \ogp+d?(M+V) logp) bit operations^ Once the norm 
has been computed the main computation is that of the p-adic logarithm, which must be done 
to the precision p M+e since we will be dividing by log p (u) G p e Z p . Since {N(f3+k-X)) G l+p e Z p 
this computation can be done using at most M multiplications of precision p M+e . Computing 
(\+X) a to the precision (p M ,X L ) takes 0~(ML logp) operations and multiplying by the inner 
product and the inverse of the norm costs 0~(MLc log p) bit operations. The result follows 
from simplifying and noting that the sum has 0(p ed d d (pM + L) d ) termsEI □ 

Corollary 5.26. With a precomputation of cost 0~(p( e+lS)d d d+3 M d+2 c) bit operations, for any 
given s in Z p we can compute 2 p (C, c; s) to the precision p M in 0~(M 2 log 2 p) bit operations. 

Proof. With L := \Mje\ we precompute 3 p (C,qX) to the precision (p M ,X L ), using the 



theorem to estimate the cost. Given s G 



we compute t 



precision p in 0~(M log p) bit operations and compute 3 P (C, c;t) 
in 0~((M 2 /e) log p) operations. The corollary follows. 



1 G p e 7L p to the 



Z P (C, c; s) (mod p 



□ 



Remark 5.27. From its definition it would seem more natural to compute 3 P (C, c; X) modulo 
(X,p e ) L . In particular, it would be enough to compute values of Z p (C, c;s). This implies 
that, for < I < L, the coefficient of X* would have to be computed to the precision p e ( L ~^ . 
By Corollary 15.151 we can replace N = p e (p\L/e] + L) with N = maxo<^<Lp e (pe(L — £) + 



20 
21 



See the proof of Theorem 15.211 for the computation of the norm. 

The costs of computing the product by A(C, 2) and the trace are negligible compared to that of the sum. 
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£) = p e+1 Le in the formula K = (N — l)d. It is clear that this does not give a significant 
improvement in the estimate of the computation time. 

We finish this subsection with a result on the direct computation of Z P (C, c; s) for a given 
s G Z p . 

Theorem 5.28. If s G Z p then 

(pM+l)d g 



Z p (C,c;s) = Tn 



A(C,E) ■^(/3 + A I -A)- s []i? feii(pM+1)(i (S(A 4 )) 



=1 



(mod p M ) 



and hence we can compute Z(C, c; s) to the precision p M in 0~(p d d d+s M d+2 c) bit operations. 

Proof. We use Lemma 15.251 again with fj, = H p q and f(x) = <f)— s (x). By Corollary 15.131 we 
can take N = pM + 2. This establishes the formula. We estimate the computation cost as in 
Theorem 15.211 replacing the computation of (1 + T) a by that of a~ s , with a := N(f3 + k ■ A), 
and accounting for the difference in the number of terms in the sum. □ 

5.4. Explicit cone decomposition. The construction of the measure fj, p 'm relies on the 
existence of a cone decomposition of a modulo m. Such a construction exists by a result 
of Cassou-Nogues [7J (see Theorem I3.6|) , based on the work of Shintani [29J , but the proof is 
nonconstructive. For d = 1 the construction is trivial, and the case d = 2 has been well studied 
(see below). For d = 3 an explicit efficient decomposition is given in [13] . but quantitative 
results on the number of discrete cones obtained at the end are missing. A general construction 
is given in [11], but this construction relies on the existence of a set of units satisfying certain 
conditions and there does not appear to be any practical way to construct such a set. However, 
the construction is generalized i n 114] for any set of units of maximal rank, at the price of using 
signed cones rather than cones This does not cause any complication in the computations 
and the changes needed to use signed cones are straightforward. 

Remark 5.29. The construction given in the present article assumes that we can always find 
a c-admissible cone decomposition. Although this is always possible in the case d = 1 and 
d = 2 (see below), it is not guaranteed by Theorem 13.61 in general. This is not really a problem 
however; one can always first construct the several cone decompositions needed, then choose 
c so that these cone decompositions are c-admissible. This is how it is done in [7J. 

We start with the case d = 1, which is straightforward. 

Proposition 5.30. Assume d = 1, that is, E = Q. L et f and a be positive integers such that 
m = /Z • oo and o = aZ and let b = a(a _1 mod /)o Then a c-admissible cone decomposition 
of a modulo m is given by the unique cone C(b;af). □ 

It is also not difficult to construct a cone decomposition in the quadratic case. Indeed, 
assume d = 2; then E is a real quadratic field. For two linearly independent elements 70 and 
71 of E + we define the half-open rational cone of a modulo m (or simply the rational cone of 
a modulo m) generated by 70 and 71 to be the following subset of a n E + : 

RC m (lo, 71; a) := { S70 + *7i with s, t e Q, < s, < t } n a n E m . 
We go from half-open rational cones to discrete cones using the formula 

■RCm(70)7i;fl) = (J C(a; 70, 71) (disjoint union) (5.9) 
w here aePC m (70,7150) 

PC m (jo, 7i; a) := { S70 + tji with s,t G Q, < s < 1, < t < 1 } D afl E m . 



That is, cones together with a ± sign; the decomposition is obtained by removing the cones with a — sign. 



2 ^In particular, by (HI), q divides /. 
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Thus from a finite family of disjoint half-open rational cones giving a set of representatives of 

afl E m modulo U m (E) we can get a c-admissible cone decomposition of o modulo m, provided 

that the generators of the rational cones are in of \ c. Let e+ be the generator of U+(E), 

the group of totally positive units of E, with > 1 > e^/. Let % be the index of U m (E) 

as a subgroup of U + (E); thus e m := e!^ is a generator of U m (E). For any totally positive 

element 7 of af \ c one can easily check that i?C m (7, e m 7; a) is a set of representatives of ofl E m 

modulo U m (E), and therefore we can construct from this rational cone a c-admissible cone 

decomposition by the formula aboveE3 However, the number of points in PC m (7, 67; 0) is of 
(2) (2) ■ 

the order of = (e + )* m and therefore much too large for computations in general. We use 
instead the following algorithm, also used in [27], which is based on the work of Hayes [19] . 
For two elements 70 and 71 of E + we set 

&(7o,7i) : = bo hi ] and #(7o,7i) := "To + 6(70,71)71 ■ 
Algorithm 5.31 (Computation of cone decomposition in degree 2). 
Input: Ideal a coprime to c and m. 

Output: A c-admissible cone decomposition of a modulo m. 

1. Compute 5 E N and h € %e such that af = Z5 + Zh. 

2. If < fcW then do h <- -h. 

Do h<r- h+ \-hW/g]g. Set (g Q , gi ) <- (g,h). 

3. While gf ] < g ( Q 2) , do (50,51) <- (51, R(go, 5i))- 

4. If g € c then do (50,51) <- (51, R(go, 5i))- 

5. Set D <- and 5i ast <- g e m . 

6. While g ^ g iast , do 

6.1. If 51 c then do 

fro <^5o, h ^-51 and (50,51) <- (51, #(5o, 5i)), 
Else 

52 «- R(go,9i), 

bo <- 5o, &i ^52 and (50,51) <- (52,^(51,52))- 

6.2. For a e PC m (b , b t ; a), do D «- D U {C(a; 6 , 61)}. 

7. Return D. 

Proposition 5.32. Lei -D^ 6e the discriminant of E and let e+ be the generator of the group 
U + (E) of totally positive units of E such that e + > 10 Then Alaorithm \5. 31\ computes a c- 
admissible cone decomposition of a modulo m in 0~ (j\f \<xf)y/D~E +A/"(f) 2 e + log De) operations 
in E. Moreover, this cone decomposition contains 0~ (N '(f) e+ log De) cones. 

Proof. The pair (5, h) constructed in Step 2 of the algorithm satisfies af = Z5 + Z/1, g^> > /i^ 1 ) 
and 1 < h^/h^ 1 '. From this pair we construct in the following steps a sequence (7 n )n>o with 
(7o,7i) := (5, h) and 7n+i := #(7n-i,7n) f° r n > 1. One can prove that the elements of this 
sequence satisfy 

(l)of = Z7 n + Z7n+i, (2)7i 1} >7Si, and (3) < T^i/t^i- 

One can also prove (see below) that there exists an integer N > 1 such that 

(4) < 7^ 2) for all n > N. 

24 We will explain below in the proof of the next proposition how to compute the elements of PC m (70, 71; o). 
2 ^To simplify the expressions, we assume from now on that E is embedded into R by the map x i-> x^ 2 \ 
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We let 7„ := ^N+n for n > 0. This is the sequence that is computed after Step 3. The 
points 7„ are successive points on the convexity polygon of of as defined in [19]. We can 
also extend the sequence in the other direction to obtain a sequence (7 n ) rag z, infinite in both 
directions, containing all the points on the convexity polygon, and for which we still have 
7 n+ i = i?(7 n _i,7 n ) for all n G Z. It will be necessary to ensure that 70 c. If 70 G c, 
we iterate one more time in Step 4 to replace 70 by 71 (that is, replacing N by N + 1). 
Indeed, by (1) 70 and 71 cannot both be in c. We can assume from now on that 70 c. 
The group U+(E) acts on the set {"fn, n £ Z}; thus there exists an integer Pq > 1 such that 
"fn+p = e +7n for all n > 0. Therefore, for any n G Z the union of the (disjoint) rational cones 
^C m (7„,7„ + i; a), . . . , RC m ("f n+ p-i,j n+ p; a), with P := i m P , gives a set of representatives 
of a n E m modulo U m (E) with generators in af. However, although 70, and thus also tp, 
do not belong to c, it is possible that ]„ £ c for some n in the range 1 < n < P — 1. 
In that case, 7„_i and 7 n +i do not lie in c, by (1), and we use in Step 6.1 the fact that 
i?C w (7 n _i,7 n ;a) U RC m (j n ,j n+1 ; 0) = i?C TO (7 n _i,7 n+1 ; 0) to get rid of cones with 7„ G c. 
We end up with a rational cone decomposition having generators suitable for constructing a 
c-admissible cone decomposition of a modulo m using (j5.9j) . 

We will now estimate the complexity of the algorithm and the number of cones obtained 
at the end. First we find an upper bound on N, the number of iterations in Step 3. Assume 
1 < t < N — 1. Since 7^ < and 7^ < 7I we have M (74+1) < M (%), and since both 
are divisible by Af(af) we have N < N (h) / W (af) + 1. Now assume, without loss of generality, 

that a/De^ = VDe > and v A Di (2) = — /De < 0. Write h = a + by/D^ with a, b G ±Z. 
Since > h^', it follows that b < 0, and also a > —b\/D E > 0, as /i G Zg. On the 
other hand, a < g — b\/DE because < g. From the fact that af = TLg + Z/i, we find that 
b = -M(af)/2g. We compute N(h) = a 2 - b 2 D E < (g - bx/D^) 2 - b 2 D E = g 2 - 2gbx/D E ~ = 
g 2 + AA(of)/Di G 0(Af(af) 2 y/D^). Thus Step 3 requires 0(J\f(af)y/D^) operations in E. 
Next we estimate the size of P = i m Po- By [191 P- 161] we have 

(2) j j (2) 

log C + log (2) £ + (1) + - log D E > P ( ° g _ £+ (1) (5.10) 
e + e + e + e + 

where 

e ( 2 ) _ e (i) A/7,y 1 



r D E nei jV(af) ' 

For n G Z we see that 6(7 n _i,7 n ) = 2 if and only if 7„ is midway between 7 ra _i and j n +\, 
in which case we say 7„ is a midpoint. If 7 n is not a midpoint we say 7„ is a vertex. If 7 n 
is a vertex we have M(^ n ) < N(o%)\/De (see [191 Prop. 5.5]) and it follows that C < eJvi ■ 
Substitution in (|5.10p gives 

P G 0(4 2) (log4 2) + log D E )) C CT^e^ log 

Lastly, we need to explain how to perform Step 6.2, we need to estimate the cost of the 
computation, and we need to estimate the number of points in the sets PC m (bo,b\; <*)■ There 
are two cases to consider, (I) (bo, b±) = (7™, 7™+i) and (II) (bo, b\) = (j n , for an arbitrary 

neZ, We have the bijection 

a/(Zb + Zbi) ^ {sb + tb x with s,t G Q, < s < 1, < t < 1} n a, 

defined as follows. For a given class a in o/(Z6o + Z&i), lift a to an arbitrary element a G a 
and write a = sbo + tb\ with s, t G Q. Then the map above sends a to 

[«Wi) : = +l)6o + (<-L*J)6i- 
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Moreover, since bo and b\ lie in af there is a well-defined map from a/(Z6o + Z&i) to a/af 
which sends a to the class of a modulo af. The set PC m (6o,6i; a) consists of precisely those 
elements [a]^^) for which a 6 a/(Z&o + Z61) is congruent to 1 modulo f. Since a and f are 
coprime this map is surjective and therefore PC m (bo, b±; a) contains exactly d/M(f) elements, 
where d := (a : Z6o + Z61). Using the methods of [9j §4.1.3], two elements «o an d ct\ of a can 
be constructed with o = Zao + 7Lu\ and a/(Z&o + Z&i) = (Z/(ioZ)d!o + (Z/diZ)ai, where cfo 
and d\ are positive integers and d := dod\. Thus the elements of a/(Z&o + Z6i) can easily be 
enumerated in 0{d) operations in E. For case (I) we have Z7 n + Z7„ + i = of and therefore 
d = M{f) and PC m (7 n , 7 n +i; a) contains only one element. For case (II) we have 

d = (a : Z7„ + Z7„ +2 ) = (a : Z7„ + Z7„ + i)(Z7 n + Z7 n+i : Z7„ + Z7„ +2 ) = Af(f) b(<y n ,-f n+1 ). 

If 7n+i is a midpoint we have 6(7™, 7n+i) = 2. We now need to estimate the size of 6(7„, Jn+i) 
when 7 n is a vertex. Since 6(7„ , 7 n +i ) = ^(7n+ p , 7n+ p +i ) it is enough to look at what happens 
for the vertices among 70, . . . , 7p -i- Writing b n := 6(7„_i,7„) to simplify the notation, we 
have by construction 

J 1 ) ~M 
7 "-i > ~(l) > llhzl 

b - 1 " b 

and therefore 

(b Po - l)(6p _! - 1) ■■■(&!- 1) lp o b Po b Po -l • • • 6l ' 



Po 



From the fact that 7p = e+70 we find that 

e ? > JJ(6, - 1). (5.11) 

The indices i for which 6j = 2, that is, corresponding to the midpoints, do not contribute to 
the product. For indices corresponding to vertices we get 



E 



1<«<P 
7i is a vertex 

We now put everything together to get the result. For Step 1, we assume that an ideal is given 
by a 2x2 integral matrix expressing a basis of the ideal of the (fixed) integral basis of E. This 
step amounts to an HNF reduction of a 2x4 matrix (see §4.7.1]). Since we can reduce the 
entries of this matrix modulo AA(af) this step takes 0~(log(A/"(af))) bit operations and hence 
is negligible. Step 3 takes 0(A/"(af)V De) operations in E. The loop in Step 6 is iterated P 
times. The most costly operation is in Step 6.2. The cost of computing ao, a±, do, and d\ 
is essentially that of an SNF reduction of a 2x2 matrix with coefficients of size < A/"(af) and 
hence can be neglected. Enumerating the elements of PC m (bo, b\; a) takes a total of 0(M(f)P) 
operations in E for the midpoints and 0(i m 7V(f)e+) for the vertices. This gives the estimate 
on the complexity of the algorithm. To count the cones we observe that the midpoints give 
O(P) cones and the vertices give 0(i m e + ) cones. The conclusion is established by using the 
fact that i m G 0(M{f)). □ 

Remark 5.33. One could ask what would happen if we were first to construct a cone decom- 
position using the algorithm without any restriction related to c, that is, deleting Step 4 and 
always doing the first part in Step 6.1, then choosing the prime ideal c so that the decom- 
position computed is c-admissible. In fact this would not change the complexity or even the 
order of the number of cones, since both are dominated by the contributions of the midpoints 
and, in fact, for these we get the same number of discrete cones in cases (I) and (II). On the 
other hand, this would probably force the norm of c to get significantly larger and that would 
adversely affect the complexity of the remaining computations. 
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Remark 5.34. The complexity and the estimate of the number of cones given by this propo- 
sition appear in practice to be very pessimistic as they are of the order of the exponential 
of Re, the regulator of E, whereas computations point towards something of the size of Re- 
Indeed, one can use f)5.11j) to show that the number of vertices among 70, . . . , 7p + is O(Re). 
However, it appears difficult to bound the number of midpoints. One can prove that if 7 n is 
a vertex then the number of midpoints following it is |_1/(1 — b n + 7^ 1 2 1 /7i 1 ^)J , and so this 
problem is related to the question of how close a quadratic irrationality can be to an integer. 
In order to bound more efficiently the number of cones one would need to bound the size of 
the b n 's. This could be done for example using (|5.1ip by finding some non-trivial lower bound 
on the number of vertices. 

5.5. Computations of p-adic L-functions. We use the results from the preceding subsec- 
tions to estimate the complexity of computing L-functions. We will make certain assumptions. 
As noted above, we assume we have computed the necessary data to work in E: ring of in- 
tegers, class group, units, etc. We assume also that we have at our disposal a prime ideal c 
satisfying the hypotheses (HI), (H2), and (H3) and the additional hypothesis 

(H4) Either x is non-trivial and x( c ) 7^ 1> ° r X 1S trivial and (c) 1 + p e+1 Z r F^ 
We assume we have computed a list of integral ideals Oj, i = 1, . . . , h m (E), coprime to c and 
m and representing all the classes of Cl m (E). Finally, we assume we have computed a cone 
decomposition for each ideal a,; we will denote by B the maximum number of cones among 
these decompositions (see the previous subsection). In what follows 5 will denote the degree 

of Q P (x)/Qp- 

Lemma 5.35. Assume the ERH. Then there exists a prime ideal c satisfying hypotheses (HI) 
through (H4), with c G 0(log 2 (AA(f)-D^)) if X is non-trivial and c G 0~(p 2m ° log 2 (-Dg;)) if x is 
trivial, where tjiq > is such that Q mo = E D Qoo • 

Proof. We use Theorem 1 of of [2]. For the case x non-trivial the application is direct. For 
the case x trivial we apply the theorem to the character p generating the group of characters 
of G&l(Ei/E). The absolute norm of the conductor of p divides the absolute norm of the 
conductor of Qm +l/Qmoj the p-adic valuation of which is v p (q) + (p m ° +1 — l)/(p — 1). The 
result follows. □ 

Theorem 5.36. Let M and N be positive integers. Under the assumptions enumerated at the 
beginning of this subsection the measures n p ]m, for i = 1, . . . , h m (E), can be computed to the 
precision (p M ,T N ) in 0~(h m (E)d d+3 BN d+1 Mclogp) bit operations. 

Proof. This follows directly from Theorem 15.211 □ 

Corollary 5.37. Let M be a positive integer and let s E Z p , with s ^ 1 if x is trivial. Under 
the assumptions enumerated at the beginning of this subsection and after precomputations of 
cost 0~(h m (E)p d+1 d d+5 BM d+2 c) bit operations, two algebraic integers /3 andj, both belonging 
to Z[x] and with ^L Pjm (x] s) lying in 1* p [x] o-nd |7L Pjm (x;s) — (3\ p < p~ M , can be computed in 
0~(h m (E)(p 2 M 3 + 5 M log p)) bit operations. Moreover, < \^\ p < 1 if x is non-trivial 

and \ j\ p = p~ e \s — l\ p if x is trivial. 

Proof. We precompute the measures //p'm, for i = 1, h m (E), to the precision (p ,T N ), 
with iV := pM + 2, the computation cost being given by Theorem 15.361 We let /3 be an 
approximation of the sum in (|4.4p to the precision p M . The values of Zp^a^ 1 , c; s), for 
i = 1, . . . , h m (E), to the precision p M are computed using Corollary 15.231 mutatis mutandis 
in 0~(h m (E)p 2 M i ) operations. The rest of the computation of (3 takes O (h m (E) 5 M log p) bit 
operations. We now let 7 be an approximation of x( c )( c ) 1_s — 1 to the precision p M . The 



'Note that we always have (c) 6 1 + p e Z p by Lemma 14.31 
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computation of 7 takes 0~(M 2 log 2 p + 5M log p) bit operations, and it follows from (|4.4p that 
7-^p,m(x; s ) nes i n %p[x] an d |7-£p,m(x> s ) ~~ P\p < £>~ M - The assertions concerning the absolute 
value of 7 are straightforward (see the proof of Corollary I4.16|) . □ 

Recall that E n Qoo = Qm and E(m) n Qoo = Q mo+mi , so that e = m + v p (q). 

Theorem 5.38. Let M and L be positive integers. Under the assumptions enumerated at the 
beginning of this subsection there exist polynomials B(X) and C{X) in Z p [x]LY], with a cost 
of 0~(h m (E)(p ed d d+3 B(pM + L) d M 2 Lc + SMlogp)) bit operations to compute, such that 



C(X)3 Pjm ( X ;X)-B(X) e { 
with ir £ O p satisfying 



p M Z p [ X ][[X}] + X L Z p [ X ][[X]] ifx is trivial or 

not of type W, 

1 -Z p [x}[[X}} + ^—Z P [ X }[[X}} otherwise, 



{X + TT * l/Vm JJ X + TT 



1 1 

7- < v n (ir) < . 

p"H-l(p-l) - PV 1 ~ p- 1 

Moreover, p-V(p-i) < |C(0)| P < 1 if X is non-trivial and C(0) = 0, |C"(0)| P = 1 ifx is trivial. 

Proof. We use the notation and results from the proof of Theorem 14, 151 The polynomial C{X) 
is an approximation modulo (p M , X L ) of the power series C(c, xi ^) an d the polynomial B(X) 
is approximation modulo (p M ,X L ) of the power series 

h m (E) 

x(a; 1 )N(a i ;X)A(a i ,c;X). 

i=l 

The first assertion follows by (|4.7p and the integral properties of 3 p , m (x; -X")) with 7r := (£— 1)/£. 
Since £ has order p m for some integer m with 1 < m < mi this proves the inequalities on 
the p-adic valuation of ir. The properties of C{X) follow from (H4). We now evaluate the 
complexity of the computation of B(X) and C(X). Let a be one of the ideals dj and let 
{Ci, . . . , Cm} be a c-admissible cone decomposition of a. Then 

m 

A(o,c;X) =^3 p (C 3 ,c;I) 
j=i 

and the computation cost of ^4(o, c; X) to the precision (p M ,X L ) follows from Theorem 15,241 
The computation time of C(X) is negligible compared to that of B{X). □ 

Corollary 5.39. Let M be a positive integer and let s E Z p; with s ^ 1 if x is trivial. Under 
the assumptions enumerated at the beginning of this subsection and after precomputations 
costing 0~(h m (E)(p( e+1 ^ d d d+3 BM d+2 Lc + 5M logp)) bit operations, two algebraic integers, j3 
and 7, both belonging to Z[%] and with 7-L p>m (x;s) lying in Z p [x] and 

[ p~ M if x is trivial or not of type W, 

can 6e computed in 0~(M 2 log p (<5/e + log p)) bit operations. Moreover, p~V(P — < |7| p < 1 
if X is non-trivial and \*y\ p = p~ e \s — l| p if x is trivial. 

Proof. We precompute the polynomials B{X) and C(X) with L := \M/e\. The precom- 
putation cost is given by Theorem 15,381 Then we compute t := n 1-5 — 1 to the precision 
p M in 0~(M 2 log 2 p) bit operations, and compute /3 to be B(t), respectively 7 to be C(t), 
to the precision p M in 0~(<5M 2 /elogp) bit operations. The result follows from (|4,6p and the 
theorem. □ 
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We conclude with the cost of computing a single value of a p-adic L-function without 
precomputations . 

Theorem 5.40. Let M be a positive integer and let s £ 7L p , with s ^ 1 if x is trivial. Under 
the assumptions enumerated at the beginning of this subsection, two algebraic integers f3 and 
7, both belonging to Z[x] and with jL p ^ m (x', s) lying in Z p [x] and \^L P)m (x] s ) ~ fi\p — P~ M > 
can be computed at a cost of 0~(h m (E)(p d d d+3 M d+2 c + 5M logp)) bit operations. Moreover, 
p-Vfa- 1 ) < |^| < i if j£ i s non-trivial and \j\ p = p~ e \s — \\ v if x ^ s trivial. 

Proof. We proceed as in the proof of Corollary 15 . 371 with the same definitions for [3 and 7. The 
cost of computing 7 is the same. We construct /3 by computing the values of Z Pjta (aJ~ , c; s), for 
i = 1, . . . , h m (E), to the precision p M , using Theorem 15.281 mutatis mutandis. The cost of this 
computation is 0~(h m (E)(p d d d+3 M d+2 c + SMlogp)) bit operations. The result follows. □ 
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